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Abstract. 

Cell motility in viscous fluids is ubiquitous and affects many biological 
processes, including reproduction, infection, and the marine life ecosystem. Here 
we review the biophysical and mechanical principles of locomotion at the small 
scales relevant to cell swimming (tens of microns and below). The focus is on 
the fundamental flow physics phenomena occurring in this inertia-less realm, 
and the emphasis is on the simple physical picture. We review the basic 
properties of flows at low Reynolds number, paying special attention to aspects 
most relevant for swimming, such as resistance matrices for solid bodies, flow 
singularities, and kinematic requirements for net translation. Then we review 
classical theoretical work on cell motility: early calculations of the speed of a 
swimmer with prescribed stroke, and the application of resistive-force theory and 
slender-body theory to flagellar locomotion. After reviewing the physical means 
by which flagella are actuated, we outline areas of active research, including 
hydrodynamic interactions, biological locomotion in complex fluids, the design 
of small-scale artificial swimmers, and the optimization of locomotion strategics. 
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1. Introduction 

Our world is filled with swimming microorganisms: The spermatozoon that fuse with 
the ovum during fertilization, the bacteria that inhabit our guts, the protozoa in our 
ponds, and the algae in the ocean. 

The reasons microorganisms move are familiar. Bacteria such as Escherichia 
coli detect gradients in nutrients and move to regions of higher concentration [1]. 
The spermatozoa of many organisms swim to the ovum, sometimes in challenging 
environments such as tidal pools in the case of sea urchins or cervical mucus in the 
case of humans [2]. Paramecium cells swim to evade predator rotifers. 

What is perhaps less familiar is the fact that the physics governing swimming at 
the micron scale is different from the physics of swimming at the macroscopic scale. 
The world of microorganisms is the world of low "Reynolds number," a world where 
inertia plays little role and viscous damping is paramount. As we describe below, 
the Reynolds number Re is defined as Re = pUL/rj, where p is the fluid density, r\ 
is the viscosity, and U and L are a characteristic velocity and length scale of the 
flow, respectively. Swimming strategies employed by larger organisms that operate at 
high Reynolds number, such as fish, birds, or insects [3, 4, 5, 6, 7, 8], do not work at 
the small scale. For example, any attempt to move by imparting momentum to the 
fluid, as is done in paddling, will be foiled by the large viscous damping. Therefore 
microorganisms have evolved propulsion strategies that successfully overcome and 
exploit drag. The aim of this review is to explain the fundamental physics upon 
which these strategies rest. 

The study of the physics of locomotion at low Reynolds number has a long history. 
In 1930, Ludwig [9] pointed out that a microorganism that waves rigid arms like oars 
is incapable of net motion. Over the years there have been many classic reviews, 
from the general perspective of animal locomotion [10], from the perspective of fluid 
dynamics at low Reynolds number [11, 12, 13, 14, 3, 15, 16], and from the perspective 
of the biophysics and biology of cell motility [17, 18, 19, 20, 21, 1]. Nevertheless, the 
number of publications in the field has grown substantially in the past few years. This 
growth has been spurred in part by new experimental techniques for studying cell 
motility. Traditionally, motile cells have been passively observed and tracked using 
light microscopy. This approach has led to crucial insights such as the nature of the 
chemotaxis strategy of E. coli [1]. Advances in visualization techniques, such as the 
fluorescent staining of flagella [22] in living, swimming bacteria, continue to elucidate 
the mechanics of motility. A powerful new contribution is the ability to measure forces 
at the scale of single organisms and single motors. For example, it is now possible 
to measure the force required to hold a swimming spermatozoon [23, 24, 25], algae 
[26] or bacterium [27] in an optical trap. Atomic force microscopy also allows direct 
measurement of the force exerted by cilia [28] . Thus the relation between force and the 
motion of the flagellum can be directly assessed. These measurements of force allow 
new approaches to biological questions, such the heterogeneity of motor behaviour 
in genetically identical bacteria. Measurements of force together with quantitative 
observation of cell motion motivate the development of detailed hydrodynamic theories 
that can constrain or rule out models of cell motion. 

The goal of this review is to describe the theoretical framework for locomotion at 
low Reynolds number. Our focus is on analytical results, but our aim is to emphasize 
physical intuition. In §2, we give some examples of how microorganisms swim. After 
a brief general review of low-Reynolds number hydrodynamics (§3), we outline the 



The hydrodynamics of swimming microorganisms 



3 



fundamental properties of locomotion without inerti (§4). We then discuss the classic 
contributions of Taylor [29], Hancock [30] and Gray [31], who all but started the 
field more than 50 years ago (§5); we also outline many of the subsequent works that 
followed. We proceed by introducing the different ways to physically actuate a flagella- 
based swimmer (§6). We then move on to introduce topics of active research. These 
areas include the role of hydrodynamic interactions, such as the interactions between 
two swimmers, or between a wall and a swimmer (§7); locomotion in non-Newtonian 
fluids such as the mucus of the female mammalian reproductive tract (§8); and the 
design of artificial swimmers and the optimization of locomotion strategies in an 
environment at low Reynolds number (§9). Our coverage of these topics is motivated 
by intellectual curiosity and the desire to understand the fundamental physics of 
swimming; the relevance of swimming in biological processes such as reproduction 
or bacterial infection; and the practical desire to build artificial swimmers, pumps, 
and transporters in microfluidic systems. 

Our review is necessarily limited to a small cross-section of current research. 
There are many closely related aspects of "life at low Reynolds number" that we do 
not address, such as nutrient uptake or quorum sensing; instead we focus on flow 
physics. Our hope is to capture some of the current excitement in this research area, 
which lies at the intersection of physics, mechanics, biology, and applied mathematics, 
and is driven by clever experiments that shed a new light on the hidden world of 
microorganisms. Given the interdisciplinary nature of the subject, we have tried to 
make the review a self-contained starting point for the interested student or scientist. 

2. Overview of mechanisms of swimming motility 

In this section we motivate our review with a short overview of mechanisms for 
swimming motility. We define a "swimmer" to be a creature or object that moves 
by changing its body shape in a periodic way. To keep the scope of the article 
manageable, we do not consider other mechanisms that could reasonably be termed 
"swimming," such as the polymerization of the actin of a host cell by pathogens of the 
genus Listeria [32], or the gas- vesicle mediated buoyancy of aquatic micoorganisms 
such as Cyanobacteria [33]. 

Many microscopic swimmers use one or more appendages for propulsion. The 
appendage could be a relatively stiff helix that is rotated by a motor embedded in the 
cell wall, as in the case of E. coli (Fig. la), or it could be a flexible filament undergoing 
whip-like motions due to the action of molecular motors distributed along the length 
of the filament, as in the sperm of many species [21] (Figs, le and f). For example, the 
organelle of motility in E. coli and Salmonella typhimurium is the bacterial flagcllum, 
consisting of a rotary motor, a helical filament, and a hook which connects the motor 
to the filament [34, 35, 36]. The filament has a diameter of w 20 nm, and traces out a 
helix with contour length ss 10 ^tm. In the absence of external forces and moments, the 
helix is left-handed with a pitch w 2.5 /im and a helical diameter w 0.5 /im [22]. There 
are usually several flagella per cell. When the motor turns counter-clockwise (when 
viewed from outside the cell body), the filaments wrap into a bundle that pushes the 
cell along at speeds of 25-35 fim/s (see §7.3.2) [37]. When one or more of the motors 
reverse, the corresponding filaments leave the bundle and undergo "polymorphic" 
transformations in which the handedness of the helix changes; these polymorphic 
transformations can change the swimming direction of the cell [22] . 

There are many variations on these basic elements among swimming bacteria. For 



The hydrodynamics of swimming microorganisms 



4 




Figure 1. Sketches of microscopic swimmers, to scale, (a) E. coli. (b) C. 
crescentus. (c) R. sphaeroides, with flagellar filament in the coiled state, (d) 
Spiroplasma, with a single kink separating regions of right-handed and left-handed 
coiling, (e) Human spermatozoon, (f ) Mouse spermatozoon, (g) Chlamydomonas. 
(h) A smallish Paramecium. 

example, Caulobacter crescentus has a single right-handed helical filament (Fig. lb), 
driven by a rotary motor that can turn in either direction. The motor preferentially 
turns clockwise, turning the filament in the sense to push the body forward [38]. 
During counterclockwise rotation the filament pulls the body instead of pushing. The 
motor of the bacterium Rhodobacter sphaeroides turns in only one direction but stops 
from time to time [39]. The flagellar filament forms a compact coil when the motor 
is stopped (Fig. lc), and extends into a helical shape when the motor turns. Several 
archaea also use rotating flagella to swim, although far less is known about the archaea 
compared to bacteria. Archaea such as the various species of Halobacterium swim more 
slowly than bacteria, with typical speeds of 2-3/xm/s [40]. Although archaeal flagella 
also have a structure comprised of motor, hook, and filament, molecular analysis of the 
constituent proteins shows that archaeal and bacterial flagella are unrelated (see [41] 
and references therein). 

There are also bacteria that swim with no external flagellar filaments. The flagella 
of spirochetes lie in the thin periplasmic space between the inner and outer cell 
membranes [42]. The flagellar motors are embedded in the cell wall at both poles 
of the elongated body of the spirochete, and the flagellar filaments emerge from the 
motor and wrap around the body. Depending on the species, there may be one or 
many filaments emerging from each end of the body. In some cases, such as the 
Lyme disease spirochete Borrelia burgdorferi, the body of the spirochete is observed 
to deform as it swims, and it is thought that the rotation of the periplasmic flagella 
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causes this deformation which in turn leads to propulsion [43, 44]. The deformation can 
be helical or planar. These bacteria swim faster in gel-like viscous environments than 
bacteria with external flagella [45, 46]. Other spirochetes, such as Treponema primitia, 
do not change shape at all as they swim, and it is thought that motility develops due to 
rotation of the outer membrane and cytoplasmic membrane in opposite senses [43, 47]. 
Finally, we mention the case of Spiroplasma, helically shaped bacteria with no flagella 
(Fig. Id). These cells swim via the propagation of pairs of kinks along the length of 
the body [48]. Instead of periplasmic flagella, the kinks are thought to be generated 
by contraction of the cytoskeleton [49, 50, 51]. 

Eukaryotic flagella and cilia are much larger than bacterial flagella, with a typical 
diameter of w 200 nm, and with an intricate internal structure [21]. The most common 
structure has nine microtubule doublets spaced around the circumference and running 
along the length of a nagellum or cilium, with two microtubules along the center. 
Molecular motors (dynein) between the doublets slide them back and forth, leading to 
bending deformations that propagate along the nagellum. There is a vast diversity in 
the beat pattern and length of eukaryotic flagella and cilia. For example, the sperm 
of many organisms consists of a head containing the genetic material propelled by a 
filament with a planar or even helical beat pattern, depending on the species [52]. 
The length of the nagellum is 12 ttm in some Hymenoptera [53], w 20 /im for hippos, 
f=s 40 /xm for humans [2] (Fig. le), s=s 80 /xm for mice (Fig. If), and can be 1mm [54] 
or even several cm long in some fruit flies [55] (although in the last case the flagella 
are rolled up into pellets and offered to the female via a "pea-shooter" effect [55]). 

Many organisms have multiple flagella. Chlamydomonas reinhardtii is an algae 
with two flagella that can exhibit both ciliary and flagellar beat patterns (Fig. lg). 
In the ciliary case, each fiagellum has an asymmetric beat pattern [21]. In the power 
stroke, each fiagellum extends and bends at the base, sweeping back like the arms of 
a person doing the breaststroke. On the recovery stroke, the fiagellum folds, leading 
as we shall see below to less drag. When exposed to bright light, the alga swims in 
reverse, with its two flagella extended and propagating bending waves away from the 
cell body as in the case of sperm cells described above [56]. Paramecium is another 
classic example of a ciliated microorganism. Its surface is covered by thousands of cilia 
that beat in a coordinated manner [57], propelling the cell at speeds of w 500/im/s 
(Fig. lh). Arrays of beating cilia arc also found lining the airway where they sweep 
mucus and foreign particles up toward the nasal passage [58]. 

3. Flows at low Reynolds number 

3.1. General properties 

We first briefly discuss the general properties of flow at low Reynolds numbers. For 
more detail we refer to the classic monographs by Happel and Brenner [59], Kim and 
Karilla [60], and Leal [61]; a nice introduction is also offered by Hinch [62], as well as 
a more formal treatment by Pozrikidis [63]. 

To solve for the force distribution on a organism, we need to solve for the flow 
field u and pressure p in the surrounding fluid. For an incompressible Newtonian fluid 
with density p and viscosity r], the flow satisfies the Navier-Stokes equations 




V-u = 0, 



(1) 
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with boundary conditions appropriate to the problem at hand. The Navier-Stokes 
equations are a pointwise statement of momentum conservation. Once u and p are 
known, the stress tensor is given by a = —pi + ^[Vu + (Vu) T ], and the force F and 
torque L acting on the body are found by integrating along its surface S 

F(t) = J J a-ndS, L(t) = J J x x (a • n) dS. (2) 

The Reynolds number is a dimensionless quantity which qualitatively captures 
the characteristics of the flow regime obtained by solving Eq. (1), and it has several 
different physical interpretations. Consider a steady flow with typical velocity U 
around a body of size L. The Reynolds number Re is classically defined as the ratio 
of the typical inertial terms in the Navier-Stokes equation, ~ pu • Vu, to the viscous 
forces per unit volume, ~ rj\7 2 u. Thus, Re = pLU/rj. A low-Reynolds number flow is 
one for which viscous forces dominate in the fluid. 

A second interpretation can be given as the ratio of time scales. The typical 
time scale for a local velocity perturbation to be transported convectively by the flow 
along the body is i a d v ~ L/U, whereas the typical time scale for this perturbation to 
diffuse away from the body due to viscosity is idiff ~ pL 2 /rj. We see therefore that 
Rc = tdiff Aadv, and a low Reynolds number flow is one for which fluid transport is 
dominated by viscous diffusion. 

We can also interpret Re as a ratio of forces on the body. A typical viscous stress 
on a bluff body is given by oviscous ~ vU/L, leading to a typical viscous force on the 
body of the form /viscous ~ vjU L. A typical inertial stress is given by a Bernoulli-like 
dynamic pressure, CTi ncrt i a i ~ pU 2 , and therefore an inertial force /inertial ~ pU 2 L 2 . 
We see that the Reynolds number is given by Re = /inertial/ /viscous j an d therefore in 
a low-Reynolds number flow the forces come primarily from viscous drag. 

A fourth interpretation, more subtle, was offered by Purcell [14]. He noted that, 
for a given fluid, T = r) 2 / p has units of force, and that any body acted upon by the 
force T will experience a Reynolds number of unity, independent of its size. Indeed, 
it is easy to see that Re = /viscous/.? 7 and Re = (/inertial/-? 7 ) 1 / 2 , and therefore a body 
with a Reynolds number of one will have /inertial = /viscous = T. A body moving 
at low Reynolds number experiences therefore forces smaller than T {T s=s 1 nN for 
water) . 

What are the Reynolds numbers for swimming microorganisms [3]? In water 
(p w 10 3 kg/m 3 , r] w 10~ 3 Pa-s), a swimming bacterium such as E. coli with 
U w 10 pm/s and L rts 1-10 /xm has a Reynolds number Re rts 10~ 5 -10~ 4 . A human 
spermatozoon with U rs 200 pm/s and L rj 50 pm moves with Re rj 10 -2 . The 
larger ciliates, such as Paramecium, have U rj 1 mm/s and L rj 100 pm, and therefore 
Re Ri 0.1 [13]. At these low Reynolds numbers, it is appropriate to study the limit 
Re = 0, for which the Navier-Stokes equations (1) simplify to the Stokes equations 

-Vp + ?,V 2 u = 0, V-u = 0. (3) 

Since swimming flows are typically unsteady, we implicitly assume the typical 
frequency u) is small enough so that the "frequency Reynolds number" pLuj 2 /rj is 
also small. Note that Eq. (3) is linear and independent of time, a fact with important 
consequences for locomotion, as we discuss below. 

Before closing this subsection, we point out an important property of Stokes 
flows called the reciprocal theorem. It is a principle of virtual work which takes a 
particularly nice form thanks to the linearity of Eq. (3). Consider a volume of fluid 
V, bounded by a surface S with outward normal n, in which you have two solutions 
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to Eq. (3), U! and u 2 . If the stress fields of the two flows are cri and er 2 , then the 
reciprocal theorem states that the mixed virtual works are equal: 



m • <T2 • n dS 1 = // U2 • <7i • ndiS. (4) 
s J Js 

3. 2. Motion of solid bodies 

When a solid body submerged in a viscous fluid is subject to a external force smaller 
than J 7 , it will move with a low Reynolds number. What is its trajectory? Since 
Eq. (3) is linear, the relation between kinetics and kinematics is linear. Specifically, 
if the solid body is subject to an external force F, and an external torque L, it will 
move with velocity U and rotation rate $7 satisfying 

F A ( A B A ( U „ (5) 



L J \ B T C I \ fl 

or the inverse relation 

U\_/M N\/F 

o r N T O ' L 



(6) 



The matrix in Eq. (5) is the "resistance" matrix of the body, and the matrix of 
Eq. (6) is the "mobility" matrix. The reciprocal theorem (4) forces these matrices to 
be symmetric [59]. Dimensionally, since low- Re stresses scale as <~ VjU/L, the sub- 
matrices scale as [A] <~ r]L, [B] ~ r]L 2 , [C] <~ 17L 3 , and similarly [M] ~ (ryL) -1 , 
[N] <~ (ryL 2 ) -1 , [O] <~ (rjL 3 )^ 1 . For most problems, the details of the geometry of the 
body make these matrices impossible to calculate analytically. The simplest example 
is that of a solid sphere of radius a, for which we have isotropic translational and 
rotational drag, M = A -1 = {(S-K-qa)^ 1 !, and O = C _1 = (87T^a 3 ) _1 l; the cross- 
couplings M and B vanish by symmetry. 

Three important features of Eqs. (5-6) needed to be emphasized for their 
implications for locomotion. The first important property is drag anisotropy: The 
matrices A, M, C, and O need not be isotropic (proportional to 1). As we discuss 
in §4, drag anisotropy is a crucial ingredient without which biological locomotion 
could not occur at low Reynolds number. For a simple illustration, consider a slender 
prolate spheroid of major axis a and minor axis b with a » b. If e denotes the 
direction along the major axis of the spheroid, we have A = ^ee + A±(l — ee), with 
A\\ w 47r?7a/[ln(2a/6) - 1/2] and A± w 8ivna[hi(2a/b) + 1/2]. 

Secondly, there exist geometries for which the matrices B and N are non-zero: 
chiral bodies, which lack a mirror symmetry plane. In that case, there is the possibility 
of driving translational motion through angular forcing — this strategy is employed by 
bacteria with rotating helical flagella (see §6). 

Thirdly, these matrices are important as they allow calculate the diffusion 
constants of solid bodies. The fluctuation-dissipation theorem states that, in thermal 
equilibrium at temperature T, the translational diffusion constant of a solid body 
is given by the Stokes- Einstein relationship D = fc^TM, where ks is Boltzmann's 
constant, while the rotational diffusion constant is given by = fc^T O. The typical 
time scale for a body to move by diffusion along its own size is td ~ L 2 /[D], while 
<~ [Dfl] -1 is the typical time scale for the reorientation of the cell by thermal forces 
alone. For a non-motile E. coli bacterium at room temperature, we have [D] ~ 0.1 
/jm 2 /s in water; while the time scale for thermal reorientation of the cell axis, [D^] _1 , 
is a few minutes. 
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3.3. Flow singularities 

Since the Stokes' equations, Eq. (3), are linear, traditional mathematical methods to 
solve for flow and pressure fields can rely on linear superposition. The Green's function 
to Stokes flow with a Dirac-delta forcing of the form S(x — x')F is given by 

1 /l rr\ 

u(x) = G(x - x') • F, with G(r) = — \- + - j , r = |r|, (7) 

p(x) = H(x - x') • F, with H(r) ^~ (8) 

The tensor G is known as the Oseen tensor, and the fundamental solution, Eq. (7), is 
termed a stokeslet [30]. Physically, it represents the flow field due to a point force, F, 
acting on the fluid at the position x' as a singularity. The velocity field is seen to decay 
in space as 1/r, a result which can also be obtained by dimensional analysis. Indeed, 
for a three-dimensional force F acting on the fluid, and by linearity of Stokes' flow, 
the flow velocity has to take the form u ~ Ff(9, r, rj), where 9 is the angle between 
the direction of F and r, and where F is the magnitude of the force. Dimensional 
analysis leads to u ~ g(6)F/r)r with a 1/r decay. 

An important property of the stokeslet solution for locomotion is directional 
anisotropy. Indeed, we see from Eq. (7) that if we evaluate the velocity in the direction 
parallel to the applied force, we obtain that = F/4-irrir, whereas the velocity in the 
direction perpendicular to the force is given by u± — F/%irrjr. For the same applied 
force, the flow field in the parallel direction is therefore twice that in the perpendicular 
direction (u|| = 2 it j_). Alternatively, to obtain the same velocity, one would need to 
apply a force in the perpendicular direction twice as large as in the parallel direction 
(F± = 2F||). Such anisotropy, which is reminiscent of the anisotropy in the mobility 
matrix for long slender bodies (§3.2; see also §5.2) is at the origin of the drag-based 
propulsion method employed by swimming microorganisms (see §4.3). 

From the fundamental solution above, Eq. (7), the complete set of singularities 
for viscous flow can be obtained by differentiation [64]. One derivative leads to 
force-dipoles, with flow fields decaying as 1/r 2 . Two derivatives leads to source- 
dipole (potential flow also known as a doublet), and force-quadrupoles, with velocity 
decaying in space as 1/r 3 . Higher-order singularities are easily obtained by subsequent 
differentiation. 

A well-chosen distribution of such singularities can then be used to solve exactly 
Stokes' equation in a variety of geometry. For example, the Stokes flow past a sphere 
is a combination of a stokeslet and a source-dipole at the center of the sphere [65]. 
For spheroids, the method was pioneered by Chwang and Wu [64], and we refer to 
Refs. [60, 66] for a textbook treatment. A linear superposition of singularities is also 
at the basis of the boundary integral method to computationally solve for Stokes flows 
using solely velocity and stress information at the boundary (see Refs. [63, 66]). 

When a flow field is described by a number of different flow singularities, the 
singularity with the slowest spatial decay is the one that dominates in the far field. 
Since a cell swimming in a viscous fluid at low- Reynolds numbers is force- and torque- 
free (Eq. 9 below), the flow singularities that describe point-forces (stokeslets) and 
point-torques (antisymmetric forcc-dipole, or rotlets) cannot be included in the far- 
field description. As a result, the flow field far from a swimming cell is in general well 
represented by a symmetric force-dipole, or stresslet [67]. Such far-field behavior has 
important consequences on cell-cell hydrodynamic interactions as detailed in §7.1. 
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4. Life at low Reynolds number 

We now consider the general problem of self-propelled motion at low Reynolds number. 
We call a body a "swimmer" if by deforming its surface it is able to sustain movement 
through fluid in the absence of external (non-hydrodynamic) forces and torques. Note 
that the "body" includes appendages such as the cilia covering a Paramecium or the 
helical flagella of E. coli. 

4-1. Reinterpreting the Reynolds number 

We first offer an alternative interpretation of the Reynolds number in the context of 
swimming motion. Let us consider a swimmer of mass m and size L swimming with 
velocity U through a viscous fluid of density p and viscosity r\. Suppose the swimmer 
suddenly stops deforming its body; it will then decelerate according to Newton's law 
ma = /drag- What is the typical length scale d over which the swimmer will coast due 
to the inertia of its movement? For motion at high Reynolds number, as in the case of 
a human doing the breaststroke, the typical drag is /drag = /inertia] ~ pU 2 L 2 , leading 
to a deceleration a ~ pU 2 L 2 /m. The swimmer coasts a length d ~ U 2 /a ~ m/(pL 2 ). 
If the swimmer has a density p s ~ m/L 3 , we see that the dimensionless coasting 
distance is given by the ratio of densities, d/L ~ p s /p. A human swimmer in water 
can cruise for a couple of meters. In contrast, for motion at low Reynolds number, 
the drag force has the viscous scaling, /drag = /viscous f]UL, and the swimmer can 
coast a distance d w LRep s /p, where p s <~ m/L 3 is the density of the swimmer. 
For a swimming bacterium such as E. coli, this argument leads to d w 0.1 nm [14]. 
For Re < 1, The Reynolds number can therefore be interpreted as a nondimensional 
cruising distance. 

A consequence of this analysis is that in a world of low-Reynolds number, the 
response of the fluid to the motion of boundaries is instantaneous. This conclusion 
was anticipated by our second interpretation of the Reynolds number (§3), where we 
saw that in the limit of very low Re, velocity perturbations diffuse rapidly relative 
to the rate at which fluid particles are carried along by the flow. To summarize, the 
rate at which the momentum of a low-Re swimmer is changing is completely negligible 
when compared to the typical magnitude of the forces from the surrounding viscous 
fluid. As a result, Newton's law becomes a simple statement of instantaneous balance 
between external and fluid forces and torques 

F ext (i)+F(i) = 0, L oxt (t) + L(t)=0. (9) 

In most cases, there is no external forces, and F ext (t) = 0. Situations where L ext (t) is 
non-zero include the locomotion of nose-heavy or bottom- heavy cells [68] ; in all other 
cases we will assume L ext (i) = 0. 

4-2. The swimming problem 

Mathematically, the swimming problem is stated as follows. Consider a body 
submerged in a viscous fluid. In a reference frame fixed with respect to some arbitrary 
reference point in its body, the swimmer deforms its surface in a prescribed time- 
varying fashion given by a velocity field on its surface, ug(i). The velocity field ug(t) 
is the "swimming gait" . A swimmer is a dcformable body by definition, but it may be 
viewed at every instant as a solid body with with unknown velocity U(i) and rotation 
rate The instantaneous velocity on the swimmer's surface is therefore given by 
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Figure 2. Physics of drag-based thrust: The drag anisotropy for slender filaments 
provides a means to generate forces perpendicular to the direction of the local 
actuation (see text for notation). 



u = U + rixx + us, which provides the boundary conditions needed to solve Eq. (3). 
The unknown values of U(i) and fl(t) are determined by satisfying Eq. (9) 

The mathematical difficulty of solving the swimming problem arises from having 
to solve for the Stokes flow with unknown boundary condition; in that regard, low- 
Re swimming is reminiscent of an eigenvalue problem. A great simplification was 
derived by Stone and Samuel [69], who applied the reciprocal theorem, Eq. (4), to 
the swimming problem. Recall that the reciprocal theorem involves two different flow 
problems for the same body. Let u and er the velocity and stress fields we seek in the 
swimming problem. For the second flow problem, suppose u and & are the velocity and 
stress fields for the dual problem of instantaneous solid body motion of the swimmer 
with velocity U and rotation rate fi. This problem correspond to subjecting the shape, 
S(t), to an external force, F, and torque, L. Applying the reciprocal theorem (4), we 
obtain [69] 

F-U + L- fi = -/7 u s ■ <5- ■ ndS. (10) 
J Js(t) 

Equation (10) shows explicitly how the swimming velocity U and rotation rate fi 
may be found in terms of the gait us, given the solution to the dual problem of the 
flow induced by the motion of the rigid body with instantaneous shape S(t), subject 
to force F and torque L. Note that since F and L are arbitrary, Eq. (10) provides 
enough equations to solve for all components of the swimming kinematics. Note also 
that for squirming motion, where the shape of the swimmer surface remains constant 
(us • n = 0), Eq. (10) simplifies further. For a spherical squirmer of radius a [69], we 
have the explicit formulas [69] 

4-3. Drag-based thrust 

Most biological swimmers exploit the motion of slender appendages ("flagella") for 
locomotion. This limit of slender bodies allows us to provide a physical, intuitive way 
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to understand the origin of locomotion through drag; the specifics of biological and 
artificial flagellar actuation will be discussed in §6. 

The fundamental property allowing for drag-based thrust of slender filaments is 
their drag anisotropy, as introduced in §3.2. Indeed, consider a thin filament immersed 
in a viscous fluid which is motionless but for flows induced by the deformations of the 
filament. The shape of the filament is described by its tangent vector t(s) at distance 
s along the filament, and its instantaneous deformation is described by the velocity 
field u(s,t), where t is time. For asymptotically slender filaments, (see §5), as in the 
case of prolate spheroids, the local viscous drag force per unit length opposing the 
motion of the filament is 

f = -£||U||-£jU±, (12) 

where U| and Uj_ are the projections of the local velocity on the directions parallel 
and perpendicular to the filament, i.e. U| = (u • t)t and ui = u — u.|| ; £|| and £j_ arc 
the corresponding drag coefficients (typically ~ 2). 

The origin of drag-based thrust relies on the following two physical ideas: (a) 
The existence of drag anisotropy means that propulsive forces can be created at a 
right angle with respect to the local direction of motion of the filament, and (b) 
a filament can deform in time-periodic way and yet create non-zero time-averaged 
propulsion. To illustrate these ideas, consider the beating filament depicted in Fig. 2. 
Any short segment of the filament may be regarded as straight and moving with 
velocity u at an angle 9 with the local tangent (Fig. 2, inset). This velocity resolves into 
components W|| = ucos9 and ui = usin#, leading to drag per unit length components 
/|| = — = — £,\\ucos6 and f±_ = — = — £±wsin#. For isotropic drag, £|| = 
and the force on the filament has the same direction as the velocity of the filament; 
however, if £|| ^ the drag per unit length on the filament includes a component 
fprop which is perpendicular to the direction of the velocity, 

Wp = (£|| ~ £_i_)wsin0cos0e x . (13) 

In addition, in order to generate a net propulsion from a time-periodic movement, 
we see from Eq. (13) that both the filament velocity u and its orientation angle 9 need 
to vary periodically in time. For example, actuation with a given u and 9, followed 
by the change u — ► — u and 9 —* n — 9, leads to a a propulsive force with a constant 
sign; in contrast, actuation in which only u — > — u changes periodically leads to zero 
average force. 

It is important to realize that this result relies on two ideas associated with 
two different length scales. One is purely local, and states that with the appropriate 
geometry and actuation, a force can be created in the direction perpendicular the 
motion of the filament. This conclusion relies explicitly on the properties of Stokes 
flows, and in a world with isotropic viscous friction (£j_ = £||), locomotion would not 
be possible [31, 70]. The second idea is a global constraint that says that the periodic 
actuation of the filament needs to be sufficiently subtle to generate non-zero forces on 
average; this property is generally known as Purcell's scallop theorem [14]. 

4-4- The scallop theorem 

As pointed out above, the Stokes equation — Eq. (3) — is linear and independent 
of time. These properties lead to kinematic reversibility, an important and well- 
known symmetry property associated with the motion of any body at zero Reynolds 
number [59, 61] . Consider the motion of a solid body with an instantaneous prescribed 
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velocity, U, and rotation rate, ft, together with the flow field surrounding it. If we 
apply the scaling U — > qU and fi — > afi, then by linearity, the entire flow and 
pressure field transform as u — > au and p — > ap. Consequently, the instantaneous 
flow streamlines remain identical, and the fluid stresses undergo a simple linear scaling, 
resulting in the symmetry F — > aF and L — > aL for the force and torque acting on 
the body. In particular, when a = — 1, this means that an instantaneous reversing of 
the forcing does not modify the flow patterns, but only the direction in which they 
are occurring. 

When applied to low-Reynolds number locomotion, the linearity and time- 
independence of Stokes equation of motion lead to two important properties [14]. 
The first one is that of rate independence: If a body undergoes surface deformation, 
the distance travelled by the swimmer between two different surface configurations 
does not depend on the rate at which the surface deformation occurs but only on its 
geometry (i.e. the sequence of shapes the swimmer is going through between these 
two configurations). 

A mathematical proof of this statement can be outlined as follows. We consider 
for simplicity swimmers with no rotational motion (an extension to the case £1 =/= 
is straightforward). Consider a body deforming its surface between two different 
configurations identified by time to and t\ . We denote by rg the positions of points 
on the surface of the swimmer. From Eq. (10), we know that the instantaneous speed 
of locomotion is given by a general integral of the form 

U(t)= // r 5 .g(r s )dS, (14) 

J J firs) 

where we have used us = rs = drs/dt. The net motion of the swimmer between to 
and ti is therefore given by 

AX = f 1 U(t)d*. (15) 

J t 

Now consider the same succession of swimmer shapes, but occurring at a different 
rate. We describe it by a mapping t' — f(t), with ti = f(U) for i = 0, 1, such that the 
shape ig(t') is the same as the shape r#(t) for all times. We now have 

AX' = r U'(t') At' = U'(/(t))/'(t) dt, (16) 

J £ J t 

where 

UW(*)= // §/'(t)-g(rs)dS= // ^-g(r s )dS(17) 

J J f{r' s ) M J J firs) M 

using the chain rule. We sec therefore that U'(t')f'(t) = U(t), and therefore 
AX = AX'. The net distance traveled by the swimmer does not depend on the rate 
at which it is being deformed, but only on the geometrical sequence of shape. One 
consequence of this property is that many aspects of low- Reynolds number locomotion 
can be addressed using a purely geometrical point of view [71, 72, 73, 74, 75]. 

The second important property of swimming without inertia is the so-called 
scallop theorem: If the sequence of shapes displayed by a swimmer deforming in a 
time-periodic fashion is identical when viewed after a time-reversal transformation 
(a class of surface deformation termed "reciprocal deformation" ) , then the swimmer 
cannot move on average. This second property puts a strong geometrical constraint 
on the type of swimming motion which will be effective at low Reynolds numbers. 
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An outline of the proof can be offered as follows (again, we consider purely 
translational motion for simplicity). Let us consider a swimmer deforming its body 
between times t and t\, and a sequence of shape described by r s . We assume that 
r s(to) = r s(ti) so that we are looking at the swimming motion over one period of 
surface deformation. The net distance traveled by the swimmer is given by Eqs. (14- 
15). We now consider the motion between ti and t% obtained by time-reversal 
symmetry of the first motion; we describe it by a temporal mapping t' = f(t), with 
ti = f(ti) and ts = f(to), defined such that the shape v' s (t') is the same as the shape 
rg(t). In that case, using similar arguments at those used to demonstrate the first 
property above, we see that 



and reversing the sequence of shape leads therefore to the opposite distance traveled. 
However, since the body deformation is reciprocal, the sequence of shape between ti 
and tz is the same as between to and t\, and therefore the distance traveled should 
be the same independently of the direction of time: AX' = AX. By combining 
the two results, we see therefore that AX = 0: Reciprocal motion cannot be used 
for locomotion at low Reynolds numbers. Note here that in order to demonstrate this 
result, we do not need to assume anything about the geometry of the fluid surrounding 
the swimmer, so the scallop theorem remains valid near solid walls, and more generally 
in confined environments. 

In his original article, Purcell illustrated this result by using the example of a 
scallop, a mollusk that opens and closes its shell in a time period fashion. A low- 
Reynolds number scallop undergoes a reciprocal deformation, and therefore cannot 
swim in the absence of inertia (independent of the rate of opening and closing)!. 
Another example of a reciprocal deformation is a dumbbell, made of two solid spheres 
separated by time-periodic distance. More generally, bodies with a single degree of 
freedom deform in a reciprocal fashion, and cannot move on average. 

Successful swimmers must display therefore non-reciprocal body kinematics. 
In his original paper, Purcell proposed a simple example of non-reciprocal body 
deformation, a two-hinged body composed of three rigid links rotating out-of-phasc 
with each other, now refereed to as Purcell's swimmer [70]. Another elementary 
example is a trimer, made of three rigid spheres whose separation distances vary in a 
time-periodic fashion with phase differences [76]. More examples are discussed in §9. 
Note that, mathematically, the presence of non-reciprocal kinematics is a necessary 
but not sufficient condition to obtain propulsion. A simple counterexample is a two 
swimmers which are mirror-images of each other and arranged head-to-hcad; although 
the kinematics of the two bodies taken together is non-reciprocal, the mirror symmetry 
forbids net motion of their center of mass. 

For biological bodies deforming in a continuous fashion, the prototypical non- 
reciprocal deformation is a wave. Consider a continuous filament of length L deforming 
with small amplitude y(x, t) (i.e. for which \dy/dx\ <C 1); in that case, the propulsive 
force generated along the filament, Eq. (13), is given by 



AX' = 




U(t) dt = - 




AX, (18) 




(19) 



| A real scallop actually swims at high Reynolds number, a regime for which the constraints of the 
theorem of course do not apply. 
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If the filament deforms as a planar wave traveling in the x- direction, y(x, t) — f(x—ct), 
the force is given by F prop — c(£|| — j f 2 (n)dr/ e x and propulsion is seen to occur 
in the direction opposite to that of the wave (oF prop < 0). Mathematically, a wave- 
like deformation allows the product (dy/dt ■ dy/dx) to keep a constant sign between 
x = and x = L, and therefore all portions of the filament contribute to generating 
propulsion. In general, all kinds of three-dimensional wave-like deformations lead to 
propulsion, in particular helical waves of flexible filaments [77]. 

Finally, it is worth emphasizing that the scallop theorem is strictly valid in the 
limit where all the relevant Reynolds numbers in the swimming problem are set to 
zero. Much recent work has been devoted to the breakdown of the theorem with 
inertia, and the transition from the Stokesian realm to the Eulerian realm is found 
to be either continuous or discontinuous depending on the spatial symmetries in the 
problem considered [78, 79, 80, 81, 82, 83, 84]. 

5. Historical studies, and further developments 

In this section we turn to the first calculations of the swimming velocities of model 
microorganisms. We consider two simple limits: (1) propulsion by small amplitude 
deformations of the surface of the swimmer, and (2) propulsion by the motion of a 
slender filament. Although these limits are highly idealized, our calculations capture 
essential physical aspects of swimming that are present in more realistic situations. 

5.1. Taylor's swimming sheet 

In 1951, G. I. Taylor asked how a microorganism could propel itself using viscous forces 
alone, rather than imparting momentum to the surrounding fluid as fish do [29]. To 
answer this question, he calculated the flow induced by propagating transverse waves 
of small amplitude on a sheet immersed in a viscous fluid. In this subsection, we 
review Taylor's calculation [29]. The sheet is analogous to the beating flagellum of a 
spermatazoon, but since the flow is two-dimensional, the problem of calculating the 
induced flow is greatly simplified. The height of the sheet over the plane y = is 

h = bsm(kx-ojt), (20) 

where the x-direction is parallel to the direction of propagation of the wave, b is the 
amplitude, k is the wavenumber, and oj is the frequency of the oscillation. Note that 
we work in the reference frame in which the material points of the sheet move up 
and down, with no x-componcnt of motion. The problem is further simplified by the 
assumption that the amplitude is small compared to the wavelength 27r/fc. Note that 
the motion of Eq. (20) implies that the sheet is extensible. If the sheet is intensible, 
then the material points of the sheet make narrow figure eights instead of moving up 
and down; nevertheless the extensible and inextensible sheets have the same swimming 
velocity to leading order in bk. 

To find the flow induced by the traveling-wave deformation, solve the Stokes 
equations with no-slip boundary conditions at the sheet, 

u(x, h(x, t)) = —bto cos(kx — Lot)e y , (21) 

with an unknown but uniform and steady flow far from the sheet, 

lim u(x,y) = -U. (22) 

y— >oo 
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Since we work in the rest frame of the sheet, U is the swimming velocity of the 
sheet in the laboratory frame, in which the fluid is at rest at y — > oo. Although it 
turns out in this problem that the leading order swimming speed is steady in time, 
other situations lead to unsteady swimming speeds. In all cases we are free to use 
non-inertial frames — even rotating frames — without introducing fictitious forces, since 
inertia may be disregarded at zero Reynolds number. 

Although U is unknown, Taylor found that no additional conditions are required 
to determine U; instead, there is a unique value of U consistent with the solution to 
the Stokes equations and the no-slip boundary condition (21). It is also important to 
note that although the Stokes equations are linear, the swimming speed U is not a 
linear function of the amplitude 6, since b enters the no-slip boundary condition both 
on the right-hand side of Eq. (21) and implicitly on the left-hand side through Eq. (20). 
In fact, symmetry implies that the swimming speed U must be an even function of b. 
Replacing b by —6 amounts to translating the wave (20) by half a wavelength. But 
any translation of the wave cannot change the swimming speed; therefore, U is even 
in b. 

Taylor solved the swimming problem by expanding the boundary condition (21) 
in bk, and solving the Stokes equations order by order. We will consider the leading 
term only, which as just argued is quadratic in b. Since the swimming velocity U 
is a vector, it must be proportional to the only other vector in the problem, the 
wavevector k = ke x . For example, if we were to consider the superposition of two 
traveling waves on the sheet, propagating in different directions, we would expect the 
swimming direction to be along the vector sum of the corresponding wavevectors. 
Dimensional analysis determines the remaining dependence of U on the parameters 
of the problem: U oc ujkb 2 . Taylor's calculation yields the proportionality constant, 
with sign: 



Note that dimensional considerations require the swimming speed to be independent 
of viscosity. This result holds due to our somewhat unrealistic assumption that the 
waveform (20) is prescribed, independent of the load. However, the rate W that the 
sheet does work on the fluid does depend on viscosity. The net force per wavelength 
exerted by the sheet on the fluid vanishes, but by integrating the local force per area 
(~ fxuibk) against the local velocity (~ cob), Taylor found 



Note that only the first-order solution for the flow is required to calculate W. In §6 we 
consider more realistic models that account for the internal mechanisms that generate 
the deformation of the swimmer. Such models can predict a viscosity-dependence 
in the swimming speed, since the shape of the beating filament may depend on 
viscosity [85, 86, 87]. And in §8 we show how the speed of a swimmer in a complex fluid 
can depend on material parameters, even for the swimming problem with prescribed 
waveform. 

According to Eq. (23), the swimmer moves in the direction opposite to the 
traveling wave. It is instructive to also consider the case of a longitudinal wave, 
in which the material points in the frame of the sheet undergo displacement S(x, t) = 
bsin(kx — ut), yielding the no-slip boundary condition 



U = -Icjkb 2 . 



(23) 



2 



W = 7juj 2 kb 2 . 



(24) 



u(x + S(x, t), y = 0) = —bui cos(kx — ut)e : 



(25) 
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Figure 3. Physical interpretation of the swimming direction for Taylor's 
swimming sheet. The vertical displacements (black vertical arrows and black 
line) along the sheet as the wave of deformation is traveling to the right (blue 
lines and rightward arrows) leads to creation of flow vorticity of alternating signs 
(curved red arrows), which induces leftward motion of the sheet (straight red 
leftward arrows). 



For a longitudinal wave, the swimming velocity is in the same direction as the traveling 
wave. 

The direction of swimming of Taylor's sheet can be understood on the basis of 
the following simple arguments. Let us consider a sheet deforming as a sine-wave 
propagating to the right (Fig. 3). Consider the vertical displacements occurring along 
the sheet as a result of the propagating wave. During a small interval of time, the 
original wave (Fig. 3, thick blue line) has moved to the right (Fig. 3, thin blue line), 
resulting in vertical displacements given by a profile which is tt/2 out of phase with 
the shape of the sheet. Indeed, wherever the sheet has a negative slope, the material 
points go up as the wave progresses to the right, whereas everywhere the sheet has 
a positive slope, the material points go down as a result of the wave. The resulting 
distribution of vertical velocity along the sheet is illustrated in Fig. 3 by the black line 
and vertical arrows. This velocity profile forces the surrounding fluid, and we see that 
the fluid acquires vorticity of alternating sign along the sheet (illustrated by the curved 
red arrows in Fig. 3). The vorticity is seen to be positive near the sheet crests, whereas 
it is negative near the sheet valleys. The longitudinal flow velocities associated with 
this vorticity distribution allow us to understand the swimming direction. In the case 
of positive vorticity, the velocity of the induced vortical flow is to the left at the wave 
crest, which is the current position of the sheet. In the case of negative vorticity, 
the induced flow velocity is to the left in the valleys, which is also where the sheet 
is currently located. As a consequence, the longitudinal flow field induced by the 
transverse motion of the sheet leads to flow contribution which are to the left in all 
cases, and the sheet is seen to swim to the left (straight red leftward arrows). Note 
that if the sheet were not free to move, then it would create an external flow field that 
cancels the sheet-induced flow, and the net flow direction would therefore be to the 
right — the sheet acts as a pump. 

There are many generalizations to Taylor's 1951 calculation. Taylor himself 
considered the more realistic geometry of an infinite cylinder with a propagating 
transverse wave [88]. In this case, there is a new length scale a, the radius of the 
cylinder, and the calculation is organized as a power series in b/a rather than bk. In 
the limit ak — > 0, the swimming velocity has the same form as the planar sheet [88]. 
With a cylinder, we can study truly three-dimensional deformations of a filament, 
such as helical waves. A helical wave can be represented by the superposition of 
two linearly polarized transverse waves, with perpendicular polarizations and a phase 
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difference of n/2. If these waves have the same amplitude, speed, and wavelength, 
then the swimming velocity is twice the velocity for a single wave. Although the 
hydrodynamic force per unit wavelength acting on the waving filament vanishes, there 
is a nonvanishing net hydrodynamic torque per unit wavelength, which is ultimately 
balanced by the counter-rotation of the head of the organism [88] . 

The Taylor sheet calculation may also be extended to finite objects. For example, 
to model the locomotion of ciliates such as Opalina and Paramecium, Lighthill 
introduced the "envelope model," in which the tips of the beating cilia that cover 
the cell body are represented by propagating surface waves [89, 90, 91]. Perhaps the 
simplest version of the envelop model is the two-dimensional problem of an undulating 
circle in the plane, which may equivalently be viewed as an infinite cylinder with 
undulations traveling along the circumferential direction [91]. Unlike the problem of a 
rigid cylinder towed through liquid at zero Reynolds number, the undulating cylinder 
does not suffer from the Stokes paradox [65, 92], since the total force on the cylinder 
is zero. And unlike the Taylor sheet problem, where the swimming speed emerges self- 
consistently, the condition of vanishing force is required to determine the swimming 
speed of the undulating cylinder. The problematic solutions that lead to the Stokes 
paradox are the same ones that lead to a net force, as well as a diverging kinetic 
energy, and are therefore eliminated in the swimming problem [91]. 

In three dimensions, the swimming speed is also determined by the condition of 
vanishing total force [89, 90] , but since the solution to the problem of towing a sphere 
with an external force is well-behaved, we may also consider solutions with nonzero 
force. These solutions must be well-behaved if we are to apply reciprocal theorem of 
§3, which gives perhaps the shortest route to calculating the swimming speed [69, 93]. 
We can also gain additional insight into why the swimming speed for a prescribed 
deformation of the surface is independent of viscosity. Using the linearity of Stokes 
flow, at any instant we may decompose the flow field generated by the swimmer into 
a "drag flow" and a "thrust flow," v = Vd + v t [3]. The drag flow is the flow induced 
by freezing the shape of the swimmer and towing it at velocity U with a force F, 
to be determined. The thrust flow is the flow induced by the swimmer's motion at 
that instant when it is prevented from moving by an anchoring force F anc h or , which 
is determined by the shape and and rate of change of shape of the swimmer at that 
instant. Superposing the two flows, and adjusting F to cancel F anc h or yields the 
swimming speed U. Since the linearity of Stokes flow implies that both F and F anc hor 
depend linearly on viscosity, the swimming velocity U does not depend on viscosity. 
Note that the same conclusion follows from an examination of the reciprocal theorem 
formula, Eq. (10). 

Finally, in the sheet calculation, it is straightforward to include the effects of 
inertia and show that if flow separation is disregarded, the swimming speed decreases 
with Reynolds number, with an asymptotic value at high Reynolds number of half the 
value of Taylor's result (23) [94, 95]. At zero Reynolds number, the effect of a nearby 
rigid wall is to increase the swimming speed as the gap between the swimmer and 
the wall decreases, for prescribed waveform [94]. However, if the swimmer operates at 
constant power, the swimming speed decreases as the gap size decreases [94]. 

5. 2. Local drag theory for slender rods 

All the calculations of the previous subsection are valid when the amplitude of the 
deflection of the swimmer is small. These calculations are valuable since they allow us 
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Figure 4. Model for a rod subject to an external force along the rod (a) or 
perpendicular to the rod (b). The arrows at the centers of the spheres represent 
external forces. The flow field of the central sphere only is shown. Each sphere 
induces a similar flow field that helps push the other spheres along. 



to identify qualitative trends in the dependence of the swimming velocity on geometric 
and, as we shall see in §8, material parameters. However, since real flagella undergo 
large-amplitude deformations, we cannot expect models based on small-amplitude 
deformations to give accurate results. Fortunately, we may develop an alternative 
approximation that is valid for large deformations by exploiting the fact that real 
flagella are long and thin. The idea is to model the flow induced by a deforming 
flagellum by replacing the flagellum with a line of singular solutions to Stokes flow 
of appropriate strength. In this subsection and the following subsection we develop 
these ideas, first in the simplest context of local drag theory, also known as resistive 
force theory, and then using the more accurate slender-body theory. 

To introduce local drag theory, we develop an intuitive model for calculating 
the resistance matrix of a straight rigid rod of length L and radius a. Our model is 
not rigorous, but it captures the physical intuition behind the more rigorous theories 
described below. Suppose the rod is subject to an external force F oxt . Suppose further 
that this force is uniformly distributed over the length of the rod with a constant force 
per unit length. Our aim is to find an approximate form for the resistance matrix, 
or equivalently, the mobility matrix, with errors controlled by the small parameter 
a/L. To this end, we replace the rod with N = L/(2a) spheres equally spaced along 
the a;-axis, with positions Xj = (2a j, 0,0). According to our assumption of uniformly 
distributed force, each sphere is subject to an external force F 0Xt /iV. If there were 
no hydrodynamic interactions among the spheres, then each sphere would move with 
velocity, u = F cx t/(67T77aiV), and the mobility matrix would be isotropic. In fact, the 
motion of each sphere induces a flow that helps move the other spheres along (Fig. 4). 

To calculate the flow u., induced by the jth sphere, recall from §3 that the flow 
induced by a moving sphere is the superposition of a stokeslet and a source-dipolc. 
Since we seek the mobility to leading order, we only keep the far-field terms; thus 

Ui(x) = — 1 (1 + e x e x ) ■ (F ext /N). (26) 

Since each sphere moves with the local flow, we identify the velocity of the ith sphere 
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with the superposition of the flow induced by the force F/N at Xj and the flows 
induced by the forces at all the other spheres: 

u(xi) = Fext/(67n7oJV) + ^u j (x i ). (27) 

Since a <C L, we replace the sum over i by an integral, noting that the density of 
spheres is l/(2a), and taking care to exclude a small segment around Xi from the 
region of integration: 

bnnaN 8iri] J \Xi — x\ N 2a 

where the prime on the integral indicates that the region of integration is from —L/2 
to L/2 except for the interval (xi — 2a, Xi + 2a). Evaluation of the integral yields 

where we have disregarded end effects by assuming \xi\ <C L. Keeping only the terms 
which are leading order in ln(L/a), and using the fact that u(xi) is constant for a 
rigid rod, we find 

u= ^ ; (l + e x e x )-f ext , (30) 

where f cx t = F cxt /(2A r a) is the externally imposed force per unit length. In our model, 
the only forces acting between any pair of spheres is the hydrodynamic force: There 
are no internal cohesive forces. Therefore, drag per unit length f = — f ext , and 

f± = -£±u±, /|| = (31) 

where _L and || denote the components perpendicular and parallel to the x-axis, 
respectively, and £j_ = 2£|| = Am] / \n(L / a) . Once again, we encounter the anisotropy 
already mentioned for slender bodies (§3.2) and stokeslets (§3.3) that is necessary for 
drag-based thrust (§4.3). 

In our derivation of Eq. (31), we assumed zero deformation since the filament 
was straight. Turning now to deformed filaments, suppose that the filament is gently 
curved, na <C 1, where k 2 = |9 2 r/<9s 2 | 2 , and r(s) gives the position of the of the 
ccntcrlinc of the filament with arclcngth coordinate s. In the limit of very small 
curvature, it is reasonable to assume that the viscous force per unit length acting on 
the curved filament is the same as the viscous force per unit length acting on a straight 
rod of the same length. Since local drag theory is an expansion in powers of 1/ ln(L/a), 
it is valid for filaments that are "exponentially thin." That is, to make l/ln(L/a) of 
order e with £ < 1, we need a/L ~ exp(— 1/e). Below in §5.3 we introduce slender- 
body theory, which has the advantage of being accurate for thin (a/L ~ e p ) rather 
than exponentially thin filaments. 

Slender-body theory also more accurately captures the hydrodynamic interactions 
between distant parts of a curved filaments. To see the limitations of the resistive force 
theory coefficients of Eq. (31), consider a rigid ring of radius R and rod diameter 2a, 
falling under the influence of gravity in a very viscous fluid. Suppose the plane of 
the ring is horizontal. Compare the sedimentation rate of the ring with that of a 
horizontal straight rod with length L = 2ttR. In both cases, each segment of the 
object generates a flow which helps push the other segments of the object down. But 
the segments of the ring are closer to each other, on average, and therefore the ring 
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falls faster. Using the coefficients of local drag theory from Eq. (31) would lead to the 
same sedimentation rate for both objects. This argument shows the limitations of our 
local drag theory. One way to improve our theory is to use a smooth distribution of 
stokeslets and source-dipoles to make a better approximation for the flow induced by 
the motion of the rod. Applying this approach to a sine wave with wavelength A leads 
to [30, 12] 

^ = uir\ * ^ 

ln(2A/a) - 5 

^ = ln(2A/a) + i- (33) 

In Ref. [12], Lighthill refined the arguments of [30] and gave more accurate values for 
and £|| . Despite the limitations of local drag theory, we will see that it is useful for 
calculating the shapes of beating flagella and the speeds of swimmers. 

In the remainder of this section, we describe some of the applications of local 
drag theory to the problem of swimmers with prescribed stroke. To keep the formulas 
compact, we work in the limit of small deflections, although local drag theory is equally 
applicable to thin filaments with large deflections. Consider the problem of a spherical 
body of radius a propelled by a beating filament with a planar sine wave [77] 

h(x,t) = bsin(kx — ut). (34) 

As in our discussion of the Taylor sheet, §5.1, we work in the frame of the swimmer. 
Thus, the problem is to find the flow velocity U that yields zero net force and moment 
on the swimmer. To simplify the discussion, we suppose that external forces and 
moments are applied to the head to keep it from rotating or moving in the y-direction. 
In real swimmers, there is a transverse component of the velocity and a rotation, which 
both play an important role in determining the swimmer's trajectory and the shape 
of the flagellum [96, 97]. 

Equation (31) gives the viscous forces per unit length acting on the filament. The 
total force per unit length has a propulsive component, Eq. (19), arising from the 
deformation of the filament, and also a drag component, arising from the resistance to 
translating the swimmer along the x direction. Integrating this force per length to find 
the total x-component of force, writing the drag force on the sphere as £,oaU = 6irr]aU, 
and balancing the force on the sphere with the force on the filament yields 

Note that only the perpendicular component of the rod velocity leads to propulsive 
thrust; the motion of the rod tangential to itself hinders swimming. Inserting the 
sinusoidal waveform (34) into Eq. (35) and averaging over a period of the oscillation 
yields 

itt\ ukb 2 

2£|| l + (£ a)/(£||.L) 

The form of Eq. (36) is similar to the result for a swimming sheet, Eq. (23); when 
L 3> a and £_l = 2£h, the two expressions are identical. For L 3> a, the swimming 
speed is independent of L for fixed k and b, since lengthening the filament increases 
the drag and propulsive forces by the same amount. 
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Figure 5. Forces on a helical segment pulled through a viscous fluid with speed 
U. Half of a helical pitch is shown. The helix is prevented from rotating by an 
external moment along x. 



Since the swimmer has finite length, we can define the efficiency as the ratio of the 
power required to drag the swimmer with a frozen shape though the liquid at speed 
U to the average rate of work done by the swimmer. In other words, the efficiency is 
the ratio of the rate of useful working to the rate of total working. To leading order 
in deflection, the efficiency e is given for arbitrary small deflection h by 

(Z\\L + Z a)U 2 

e = — — j—. (37) 

Uf (h 2 )dx 

For the sinusoidal traveling wave (34) with £j_ = 2£i i , 
_ k 2 b 2 /2 

e i + ( } 

For small deflections, kb <C 1, and the hydrodynamic efficiency is small. Note that 
the total efficiency is given by the product of the hydrodynamic efficiency and the 
efficiency of the means of energy transduction. 

We now turn to the case of rotating helix, such as the flagellar filament of E. 
coli. The body of the cell is taken to be a sphere of radius a. For simplicity, suppose 
the radius R of the helix is much smaller than the pitch of the helix, or equivalently, 
that the pitch angle a is very small (Fig. 5). Expressions relevant for large amplitudes 
may be found in Rcfs. [98, 27]. We also assume that the axis of the helix always 
lies along the a:-axis, held by external moments along y and z if necessary. The 
helix is driven by a rotary motor embedded in the wall of the body, turning with 
angular speed f2 m relative to the body. The helix rotates with angular speed f2 in 
the laboratory frame, and the body must counter-rotate with speed Qb to ensure the 
total component of torque along x of the swimmer vanishes. The angular speeds are 
related by fib + fi = fi m - 

To calculate the force and moment acting on the helical filament, use Eq. (5), 
simplified to include only components along x-axis: 



F \ A B\ U 

M [ B c \ fi 



(39) 



where F and M are the external force and moment, respectively, required to pull the 
helix with speed U and angular rotation rate fi. To leading order in a, and with the 
assumption that L> a,R, the resistance coefficients are approximately 

A^^L (40) 

B~ -(Z±-Z\\)aRL (41) 

C w UR 2 L- (42) 
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These values are deduced by applying Eqs. (31) to special cases of motion. The 
coefficient A may be found by considering the case of a helix pulled at speed U but 
prevented from rotating. Likewise, C may be found by examining a helix that rotates 
about x but does not translate along x. To find B, we may examine the moment 
required to keep a helix pulled with speed U from rotating (Fig. 5). Or we may also 
examine the force required to keep a helix rotating at angular speed f2 from translating. 
The equivalence of these two calculations is reflected in the symmetry of the resistance 
matrix, which ultimately stems from the reciprocal theorem, Eq. (4). Note that the 
sign of the coupling B between rotational and translational motion is determined by 
the handedness of the helix. 

To find the swimming speed U and rate of filament rotation, we equate the 
external forces and moments on the filament to the forces and moments acting on 
the cell body, 



where we have introduced a resistance coefficient £ r = 87r for the rotation of the sphere. 
Solving for the three unknowns U, 0, and fib, we find 



fib = Q m + (D(a/L), and fi w fi m (£r/£±)a 3 /(R 2 L). The resistance £o of the body does 
not enter since we assume L»n. The speed is linear in a since the sign of the speed 
is given by the handedness of the helix, which is given by the sign of a. Note the 
contrast with the planar wave. If a = 0, then the velocity (45) vanishes: Propulsion 
by means of a rotating helix requires a body. However, a planar wave does not require 
a body for propulsion [see Eq. (36)]. For the helical filament, the swimming speed 
decreases with increasing L for fixed body size a and fixed motor speed fi m , since 
torque balance forces the filament to rotate more slowly with increasing L. Also, since 
the resistance to rotation of a helix scales as R 2 , the helix rotates sufficiently faster 
as the radius decreases to overcome the reduced thrust force implied by the linear 
dependence of B on R, leading to a speed that increases like 1/R for decreasing R. It 
is important to note that the motor speed fi m of a swimming bacterium is not directly 
observable with current techniques; typically the approach just described is used to 
make a prediction for the relation between observables such as U and fi, which is then 
compared with measurements [27, 99] 

5.3. Slender-body theory 

The local drag theory illustrated in the previous section allows an intuitive 
presentation of the scaling laws for flagella-based locomotion. It turns out however 
that this theory is quantitatively correct only for exponentially slender filaments. Let 
A be the typical length scale along the flagcllum on which its variations in shape 
occur, such as the wavelength, and a the flagellum radius. The local drag theory 
assumes that l/log(A/a) <C 1. Since real biological flagclla have aspect ratios on 
the order a/A <~ 10~ 2 , an improved modeling approach is necessary [12, 100]. The 
idea, termed slender-body theory and pioneered by Hancock [30] , is to take advantage 
of the slenderness of the filaments and replace the solution for the dynamics of the 
three-dimensional of the filament surface by that of its centerline using an appropriate 



F = - hall 
M = - £ r a 3 fi b , 



(43) 
(44) 




(45) 
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distribution of flow singularities. Two different approaches to the method have been 
proposed. 

The first approach consists of solving for the flow as a natural extension of the 
local theory, and approximating the full solution as a series of logarithmically small 
terms [101, 102, 103]. Physically, the flow field close to the filament is locally two- 
dimensional, and therefore diverges logarithmically away from the filament because 
of Stokes' paradox of two-dimensional flows [61]. The flow far from the filament is 
represented by a line distribution of stokcslets of unknown strengths, and diverges 
logarithmically near the filament due to the singular nature of the stokeslets (the 
logarithmic divergence is due to the line integration of \ jr stokeslet terms). Matching 
these two diverging asymptotic behaviors allows the determination of the stokeslet 
strengths order by order as a series of terms of order 1/ (log A/a)™. The leading-order 
term in this series, of order 1/ log(A/a), is the local drag theory, and gives a stokeslet 
distribution proportional to the local velocity. The next order term is in general non- 
local and provides the stokeslet strength at order 1/ (log A/a) 2 as an integral equation 
on the filament shape and velocity. Terms at higher order can be generated in a 
systematic fashion [101, 102, 103]. This approach to slender-body hydrodynamics is 
the logical extension to the local drag theory, and all the terms in the expansion can 
be obtained analytically which makes it appealing. However, the major drawback to 
this approach is that each term in the expansion is only smaller than the previous term 
by a factor 1/ (log A/a), so a large number of terms is necessary in order to provide an 
accurate model for the flow. 

A second approach, asymptotically more accurate but technically more involved, 
consists in bypassing the logarithmically-converging series by deriving directly the 
integral equation satisfied by the (unknown) distribution of singularities along 
the filament. Such approach has been successfully implemented using matched 
asymptotics [104] or uniform expansions [105], and leads to results accurate at order 
a I A. This approach is usually preferable since instead of being only logarithmically 
correct it is algebraically correct. This improved accuracy comes at a price however, 
and at each instant an integral equation must be solved to compute the force 
distribution along the filament. An improvement of the method was later proposed 
by accurately taking into account end effects and a prolate spheroidal cross-section, 
with an accuracy of order (a/A) 2 log(a/A) [106]. 

In his John von Neumann lecture, Lighthill proposed an alternate method for 
the derivation of such integral equations. Instead of using asymptotic expansions, he 
used physical arguments to derive the type and strength of the singularities located 
along the filament [12, 107]. By analogy with the flow past a sphere, Lighthill first 
proposed that a line distribution of stokeslets and source-dipoles should be appropriate 
to represent the flow field induced by the motion of the filament. He then demonstrated 
that the strength of the dipole distribution should be proportional to the stokcslets 
strengths using the following argument (see Fig. 6). Consider a location sq along the 
filament. By assumption of slenderness, it is possible to find an intermediate length 
scale q along the filament such that a <C q <C A. The flow field on the surface of 
the filament at the position so is then the sum of the flow due to the singularities 
within a distance q from sq ("inner" problem) and those further away than q from 
so ("outer" problem). Since ij»a, the contribution at so from the outer problem 
is primarily given by a line distribution of stokeslets (the source dipoles decay much 
faster in space). In the inner problem, Lighthill then showed that it was possible to 
analytically determine the strength of the dipoles to ensure that the complete solution 
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Figure 6. LighthilFs singularity construction for the slender-body theory of 
flagellar hydrodynamics. To compute the flow field at a point so along the 
flagellum, Lighthill picks an intermediate length scale q such that a « 9 C J 
and represents the flow by a distribution of stokeslests at a distance further away 
than q from so (outer problem) and a distribution of stokeslets and source-dipoles 
within q of so (inner problem). The dipole strengths in the inner problem are 
obtained by requiring that the result be independent of q. 



(sum of the inner and outer problems) was independent of the value of q. The dipole 
strength is found to be proportional to the stokeslet strength, and the resulting value 
of the velocity of the filament at s = s is given by the integral equation 

u(so ) = + f G(r -r)-f( S )d S . (46) 

In Eq. (46), f is the local strength of the (unknown) stokeslet distribution (dimensions 
of force per unit length), G is the Oseen tensor from Eq. (7), 5 is a length scale that 
appears in the analysis (<5 = ^ay/e), and fj_ represents the normal component of the 
stokeslet distribution, i.e. f± = (1 — tt) ■ f if t is local the tangent to the filament. 
Note that Lighthill's slender-body analysis is less mathematically rigorous than those 
presented in Refs. [104, 105], and consequently gives results which are only valid at 
order (a/ A) 1 / 2 [3]. His derivation provides however important physical insight into a 
subject, the topic of flow singularities, that is usually very mathematical. The resulting 
integral formulation, Eq. (46), is relatively simple to implement numerically, and can 
also be used to derive "optimal" resistance coefficients for the local drag theory (see 
§5.2). His modeling approach has been extended for filament motion near a solid 
boundary [108, 109, 110, 111], and an alternative approach based on the method of 
regularized flow singularities [112] has also been devised [113]. 

6. Physical actuation 

In this section we describe common mechanisms for the generation of non-reciprocal 
swimming strokes. In addition to biological mechanisms such as rotating helices or 
beating flagella, we also describe simple mechanisms that do not seem to be used by 
any organism. The non-biological setups are useful to study since they deepen our 
understanding of the biological mechanisms. For example, the modes of an elastic rod 
driven by transverse oscillations at one end are useful for understanding the shape of a 
beating flagellum driven by motors distributed along its entire length. An important 
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Figure 7. The torque-speed characteristic for the flagellar motor of E. coli for 
three different temperatures (bold lines). Also shown (thin lines) are the load 
lines for a cell body with radius of 1 fim (left), and a miniccll with effective radius 
of about 0.3 (im (right). Figure reprinted from Ref. [35], courtesy of Howard Berg 
and by permission from Annual Reviews, copyright 2003. 



theme of this section is the fluid-structure interaction for thin filaments in viscous 
liquid. 

6.1. Boundary actuation 

We begin with the case of boundary actuation, in which an elastic filament is driven by 
a motor at its base. The rotary motor of the bacterial flagellum is a prime example of 
such a biological actuating device [34, 35, 114, 115]. The steady-state relation between 
motor torque and motor speed is shown in Fig. 7 [35] . At low speeds, the motor torque 
is roughly constant; at higher speeds it decreases linearly with speed, reaching zero 
torque at about 300 Hz at 23° C. To determine the speed of the motor from the motor 
torque-speed relation, use torque balance and equate the motor torque with the load 
torque. By the linearity of Stokes flow, the load torque is linear in rotation speed. In 
the experiments used to make the graph of Fig. 7, the flagellar filament of E. coli was 
tethered to a slide, and the rotation of the body was observed. A typical body of 1 /im 
radius has a substantial resistance, leading to the steep load curve on the left of Fig. 7 
and a correspondingly low rotation speed. A smaller load, such as that of a minicell, 
leads to a load curve with smaller slope, and higher rotation speed. The torque- 
speed characteristic M(fl m ) allows us to go beyond the artifice of the previous section 
where we calculated the swimming speed U in terms of the motor speed f2 m . Solving 
Eqs. (39) and (43-44) along with M — M(fi m ) yields the swimming speed in terms of 
the geometrical parameters of the flagellum, the cell body, the drag coefficients, and 
the properties of the motor. 

In the previous section we described how a rotating helix generates propulsion. 
Since the flagellar filaments of E. coli and S. typhimurium are relatively stiff, a helical 
shape is necessary to escape the constraints of the scallop theorem as a straight 
rod rotating about its axis generates no propulsion. Indeed, mutant E. coli with 
straight flagella do not swim [116]. If the rate of rotation of a straight but flexible 
rod is high enough for the hydrodynamic torque to twist the rod through about 
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Figure 8. A bent elastic rod. The darkly shaded part of the rod exerts a force 
F(a;) and moment M(x) on the lightly shaded part through the cross-section at 
x. 



one turn, then the rod will buckle into a gently helical shape that can generate 
thrust [117, 118, 119]. However, the high twist modulus of the filament and the 
low rotation rate of the motor make this kind of instability unlikely in the mutant 
strains with straight flagella. On the other hand, a rotating helix with the dimensions 
of a flagellar filament experiences much greater hydrodynamic torque since the helical 
radius (microns) is much greater than the filament radius (w 10 nm). The handedness 
of the helix also breaks the symmetry of response to the sense of rotation of the motor: 
Counterclockwise rotation of a left-handed helix in a viscous fluid tends to decrease the 
pitch of a helix, whereas clockwise rotation tends to increase the pitch [120]. There is 
no noticeable difference between the axial length of rotating and de-energized flagella 
for counterclockwise rotation [121]; calculations of the axial extension [120] based on 
estimates of the bending stiffness of the flagellar filament [122] are consistent with this 
observation. However, the hydrodynamic torque for clockwise rotation is sufficient to 
trigger polymorphic transformations, in which a right-handed helical state invades the 
left-handed state by the propagation of a front [123, 124]. 

Now consider the case of an elastic filament driven by a mechanism that oscillates 
the base of the filament in the direction normal to the tangent vector of the filament. 
Although we know of no organism that uses this mechanism to swim, study of 
this example has proven instructive. In early work, Machin [125] pointed out that 
the overdamped nature of low Re flow leads to propagating waves of bending with 
exponential decay of the amplitude along the length of the filament. Since the 
observed beating patterns of sperm flagella typically have an amplitude that increases 
with distance from the head, Machin concluded that there must be internal motors 
distributed along the length of the flagellum that give rise to the observed shape. 
This problem has served as the basis for many subsequent investigations of the 
fluid-structure interaction in swimming [125, 126, 127, 128, 129, 130, 97, 131, 87], 
and has even been applied to the determination of the persistence length of actin 
filaments [132, 133]. Therefore, we give a brief overview of the most important elements 
of the problem. 

Consider a thin rod of length L constrained to lie in the xy-pl&ne, aligned along 
the cc-axis in the absence of external loads. We will consider small deflections h(x) 
from the straight state. When the rod is bent into a curved shape, the part of the rod 
on the outside of the curve is under tension, while the part of the rod on the inside of 
the curve is under compression (Fig. 8). Therefore, the section of the rod to to right 
of x exerts a moment M = Ak(x) w Ad 2 h/dx 2 on the section to the left of x, where 
A is the bending modulus and k(x) is the curvature of the rod at x [134]. Working 
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to first order in deflection, balance of moments on an element of length dx of the rod 
implies 

where F y = —Ad 3 h/dx 3 is the y-componcnt of the force exerted through the cross- 
section at x. Thus, if the rod has a deflection h(x), then an elastic force f y dx = 
(dF y /dx)dx = — (Ad 4 h/dx 4 )dx acts on the element of length dx at x. Balancing 
this elastic force with the transverse viscous force from resistive force theory yields a 
hyper-diffusion equation 

, dh A d 4 h 

^1h = - a W (48) 

The shape of the rod is determined by solving Eq. (48) subject to the appropriate 
boundary conditions, which are typically zero force and moment at the far end, x = L. 
At the near end, common choices are oscillatory displacement h(0,t) = bcos(u>t) with 
clamping dh/dx\ x=0 = 0, or oscillatory angle dh/dx\ x=a = 9cos(ojt) with h(0,t = 0), 
where b/L < 1 and 9 < 1 [125, 127]. 

The appearance of dh/dt in Eq. (48) causes the breakdown of kinematic 
reversibility: Even for a reciprocal actuation such as h(0,t) = bcos(uit), Eq. (48) 
implies that the rod shape is given by propagating waves. Physically, the breakdown 
of kinematic reversibility occurs because flexibility causes distant parts of the rod to 
lag the motion of the rod at the base. We saw in §3 that zero-Re flow is effectively 
quasistatic since the diffusion of velocity perturbations is instantaneous. When the 
filament is flexible, the time it takes for perturbations in shape to spread along the rod 
scales as £±L 4 /A. Since the shape of the rod does not satisfy kinematic reversibility, 
the flow it induces does not either, and a swimmer could therefore use a waving elastic 
rod to make net progress. 

The wavelength and the decay length of the propagating waves is governed 
by a penetration length, I = [A/ (^^oj)] 1 / 4 . Sometimes this length is given in the 
dimensionless form of the "Sperm number," Sp = L/C = L(cj^/^4) 1 / 4 . If the rod is 
waved rapidly, the penetration length is small I <C L and propulsion is inefficient since 
most of the filament has small defection and contributes drag but no thrust. For small 
frequencies, the rod is effectively rigid L, and there is no motion since kinematic 
reversibility is restored. Thus, we expect the optimum length for propulsion is I s=a L, 
since at that length much of the rod can generate thrust to compensate the drag of 
pulling the filament along x [127, 97]. Note that our discussion of flexibility may be 
generalized to other situations; for example, the deformation of a flexible wall near a 
swimmer is not reversible, leading to a breakdown of the scallop theorem even for a 
swimmer that has a reciprocal stroke [135]. In this case the average swimming velocity 
decays with a power of the distance from the wall, and therefore this effect is relevant 
in confined geometries. 

An interesting variation on the Machin problem is to rotate a rod which is tilted 
relative to its rotation axis. If the rod is rigid, then it traces out the surface of a 
cone. But if it is flexible, then the far end will lag the base, and the rod has a helical 
shape. As long as the tilt angle is not too small, this shape may be determined 
without considering effects of twist [136, 137]. If the driving torque rather than speed 
is prescribed, there is a transition at a critical torque at which the shape of the rod 
abruptly changes from gently helical to a shape which is much more tightly wound 
around axis of rotation, with a corresponding increase in thrust force [138, 139, 140]. 
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Figure 9. Cross-section section of an axoneme from wild type Chlamydomonas, 
courtesy of Steve King. The diameter is approximately 200 nm. Note that the 
plane containing the active motors is perpendicular to the plane defined by the 
beating filament. 



6.2. Distributed actuation 

We now consider distributed actuation, in which molecular motors arc distributed 
along the length of the filament. Eukaryotic flagella and cilia use distributed actuation. 
Figure 9 shows a cross-section of the axoneme, the core of a eukaryotic flagcllum. As 
mentioned in the introduction, the axoneme consists of nine microtubule doublets 
spaced along the circumference of the flagellum, with two microtubules running along 
the center. In this review we restrict our attention to the case of planar beating, 
although many sperm flagella exhibit helical beat patterns, and nodal cilia have a 
twirling, rotational beat pattern [141, 142]. The bending of the eukaryotic flagellum 
arises from the relative sliding of neighboring microtubule doublets [143, 144, 145, 146]. 
The sliding is caused by the action of ATP-driven dynein motors, which are spaced 
every 24 nm along the microtubles [85] . Since the relative sliding of the microtubules at 
the end near the head is restricted [147], and since each microtubule doublet maintains 
its approximate radial position due to proteins in the core of the flagellum, the filament 
must bend when the motors slide microtubule doublets. For example, in Fig. 10, 
motors have slid the lower doublet to the right relative to the upper doublet. 

The simplest approach to understanding how the sliding of the microtubules 
generates propulsion is to prescribe a density of sliding force and deduce the 
shape of the flagcllum and therefore the swimming velocity from force and moment 
balance. This approach is taken in Refs. [148, 87], where the effects of viscosity 
and viscoelasticity are studied. A more complete model would account for how 
the coordination of the dynein motors arises. Over the years, several different 
models for this coordination have been suggested. Since sea urchin sperm flagella 
continue to beat when they have been stripped of their membranes with detergent 
[149], it is thought that the motor activity is not coordinated by a chemical 
signal but instead arises spontaneously via the mechanics of the motors and their 
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Figure 10. Simple model for a planar flagellum actuated by sliding filaments. 
The dotted line is the center of the flagellum when it is straight. The solid lines 
represent the two microtubule doublets that the active motors slide relative to 
each other. The two solid curves and the dotted line all have the same length; 
the dots divide each curve into quarters as measured along the contour. 



interaction [150, 151]. A detailed discussion of molecular motors and the different 
mechanisms that have been put forth for controlling the beat pattern would take us 
too far afield [152, 153, 154, 155, 156]. Instead we review regulation by load-dependent 
motor detachment rate as presented by Riedel-Kruse and collaborators, who showed 
that this mechanism is consistent with observations of the flagellar shape [147]. 

The first step is to simplify the problem by considering planar beating and 
projecting the three-dimensional flagellum of Fig. 9 onto a planar flagellum such as 
that shown in Fig. 10. Since the lateral spacing between the two filaments is a constant, 
a, the amount that the top filament slides by the bottom is given by 

A(x) = A + a [0(x) - 9(0)} , (49) 

where Ao is the sliding displacement at the base of the flagellum where x — 0, and 9(x) 
is the angle the tangent of the flagellum centerline makes with the ir-axis. In Fig. 10, 
Ao = 0, and A < except at either end of the flagellum. We will forbid sliding at the 
base, but it is straightforward to allow such sliding in our model, and measurements 
of the shape of bull sperm flagella suggest that A(0) ~ 55 nm [147]. Equation (49) 
in the case Ao = is most easily deduced when the two filaments form the arcs of 
concentric circles, but it holds for more general curved shapes as well. Also, we will 
work in the limit of small deflection, but Eq. (49) is valid even for large deflections. 

Now consider the forces between the two filaments. Let f(x,i) denote the force 
per unit length along x that the bottom filament exerts on the top filament. This 
force per unit length could arise from passive resistance to sliding as well as motors. 
The internal moment acting a cross section has a contribution from bending and an 
additional piece from /: 

8 2 h f L 

M=A dtf- a J a fdX - (50) 
Moment balance on an element of the flagellum yields the transverse force acting on 
a cross section, F y — —Ad 3 h/dx — af, which with force balance yields the equation 
of motion for the filament: 

^i-^-f , , <"> 

For example, if the motors are not active and there is only elastic resistance to sliding, 
/ = —KA, then the shape equation is 

& dh A 8 4 h „ 2 8 2 h 

^M=- A lte +Ka l^- (52) 
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Even if the bending stiffness A vanishes, the elastic resistance to sliding will damp 
out any perturbation of the straight shape. More generally, the passive resistance will 
have a viscous as well as an elastic component: / = —KA — AA. 

Now consider an active component to /. As in the case of the bacterial rotary 
motor, the speed of a dynein motor is thought to decrease with load. One end of the 
dyncin is strongly attached to a microtubule doublet, while the other end attaches to 
and detaches from the neighboring doublet. The motor only does work when attached. 
The proposal of Refs. [151, 85, 86] is that oscillations arise spontaneously since the 
detachment rate increases with load. To see that a load-dependent detachment rate 
leads to positive feedback, consider a collection of motors sliding one filament past 
another. If there is a perturbation in which the sliding rate increases, then the 
load on each individual motor must decrease, according to the force-speed motor 
characteristic. On the other hand, the decrease in force per motor leads to an 
reduction in the detachment rate, which in turn leads to a greater total sliding force, 
and ultimately, an increase in sliding velocity. Working in the frequency domain, we 
may define a susceptibility \ y i a / — ~x( w )A. For the passive elements mentioned 
above, x = K + A simple quantitative model with load-dependent detachment 
rate leads to a susceptibility \ that can have a negative effective elastic constant or 
a negative effective drag term [147]. By using Eq. (51) to relate the shape of the 
beating flagellum to the susceptibility \, Riedel-Kruse and collaborators showed that 
the passive components of \ were small compared to the active components; in other 
words, the forces and moments induced by the motors are balanced by viscous drag and 
bending forces, rather than internal resistance to sliding. Since the oscillations arise 
spontaneously in this model, the calculation of the beating filament shape amounts to 
solving an eigenvalue problem [85, 86]. 

We close this section by noting that there are several other distinct means of 
actuation. The bacterium Spiroplasma has an internal helical ribbon that is though 
to undergo contractions that in turn cause the body to change shape. The shape 
change amounts to the propagation of kinks down the body, where the kinks are 
generated by a change in the handedness of the cell body [48, 157]. Spirochetes 
such as Leptospiracaeae have internal flagella that wrap around the periplasmic space 
between the cell body and an outer sheath. The flagellar filaments emerge from motors 
at either pole of the cell, but the rotating flagella distort the body in a nonreciprocal 
way that leads to locomotion [158, 159]. These swimmers have inspired recent table- 
top experiments in which rigid supcrhcliccs sediment in a very viscous fluid [160]. The 
helical shape causes the superhelices to rotate as they fall. An important finding of 
Ref. [160] is that resistive force theory give poor quantitative results for the rotation 
speed, even getting the sense of rotation wrong for tightly coiled superhelices. 

7. Hydrodynamic interactions 

7.1. Interactions between cells 

Microorganisms swimming in viscous fluids typically do so within semi-dilute or 
dense cell populations. As an organism is swimming, it sets up a flow which will 
be felt by the cells nearby, possibly affecting the dynamics at the level of the 
entire population. For example, spermatozoa involved in human reproduction may 
swim in population sizes as high as millions of cells [2]. Bacterial suspensions 
are known to display so-called "bacterial turbulence," where large-scale intermittent 
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Figure 11. The flow field created by a swimmer at low Reynolds number, (a): 
Cells which are pushers have a positive force-dipole (p > 0, see text) and induce a 
flow field directed away from the cells along their swimming direction (repulsion) 
and a flow field directed toward the cells along their side (attraction) (red solid 
arrows represent local forcing from the cell on the surrounding fluid); (b): Pullers 
have a negative force dipole (p < 0), inducing an attractive flow field along their 
swimming direction and a repulsive flow field along their side; (c) Two pushers on 
a converging course reorient each other, tending toward a configuration with cells 
parallel and swimming side-by-side (h is the separation distance between the cells, 
and 6 is the angle between the cell direction and the direction normal to their 
separation); (d) Two pullers on a diverging course reorient each other, tending 
toward a configuration in which the cells are antiparallel, swimming away from 
each other. 



motion in the forms of swirls and jets is set up when the cells become sufficiently 
concentrated [161, 162, 163, 164, 165, 166, 167]. Even for small numbers of cells, 
hydrodynamic interactions are suspected to play an important role, in particular in 
reproduction. Such is the case for spermatozoa of the wood mouse which aggregate, 
and thereby swim faster [168]. The pairing of the opossum spermatozoa enable them 
to swim more efficiently in very viscous fluids [169]. Fishfly spermatozoa cluster in 
dense bundles for similar reasons [170]. Recently, sea urchin sperm cells were observed 
to arrange into periodic vortices [171]. 

Let us first consider the interactions between two swimming cells. The flow field 
created by the first organism (cell A) will in general lead to two types of passive 
hydrodynamic interactions with its neighbor (cell B). First, cell B will feel the velocity 
field created by cell A, and will be carried along by this flow as a result. In addition, 
cell B will also feel the gradients in the velocity field created by cell A, which tend to 
change the orientation of B and and thereby affect its future swimming direction. 

The first type of hydrodynamic interaction can be intuitively understood by 
considering the far-field flow created by a swimming cell (Fig. 11). As discussed 
in §3.3, since a cell is force-free, the velocity in the far field is a force dipole decaying 
typically as 1/r 2 . Two different types of force-dipoles can in general arise, leading to 
significantly different physical pictures. Consider a microorganism with an elongated 
body and let e be a unit vector attached to the cell and pointing along the swimming 
direction, which is also along the elongated direction of the cell. Usually, the force 
dipole p will also be in the swimming direction, p = pe, but it can have either a 
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positive or negative sign. Cells for which p > are called "pushers," and include 
typical swimming spermatozoa, or bacteria such as E. coli. Pushers repel fluid from 
the body along the long axis, and draw fluid in to the sides (Fig. 11a; the red solid 
arrows represent local forcing from the cell on the surrounding fluid) . Cells with p < 
are called "pullers," the prototypical example of which is the alga Chlamydomonas (see 
also the artificial swimmer in Ref. [172]). Pullers draw fluid in along the elongated 
direction, and push fluid out from the sides (Fig. lib). Mathematically, the flow 
induced at a distance r from the dipole p — pe is given by 

U W = ^3[ 3COs2 ^- 1 ] r ' 

where 9 is the angle between the dipole direction, e, and the flow position, r. 
Physically, the dipole strength scales as \p\ ~ r/UL 2 where U is the swimming 
speed of the cell, and L its typical length. From Eq. (53), we see that two similar 
cells swimming side by side (9 = n/2) experience a relative velocity scaling as 
Aw s idc ~ —p/r/r 2 . Therefore two pushers attract each other, while two pullers repel 
each other [173]. Conversely, the relative velocity of the two swimmers are aligned 
along their swimming direction (9 = 0) scales as Au a ii gnc( j ~ p/^r 2 , and the sign of 
the resulting attraction/repulsion is reversed. 

The second effect of hydrodynamic interactions between two cells is the 
reorientation due to velocity gradients. The vorticity field induced by the flow set 
up by the cell, Eq. (53), is given by 

_ 3p (e • r)(e x r) 
u = Vxu=- J -- J -\ -• 54) 

Airn r° 

A sphere subject to this flow will rotate at a rate of half the vorticity, fi = ^w, to 
leading order in the ratio of the sphere size and the separation of the sphere from the 
microorganisms inducing the flow [60]. For interacting cells, which are typically not 
spherical, an additional component for f2 arises from the symmetric part of the rate of 
strain, E = ^(Vu + Vu T ), since elongated cells tend to align with the principal axis 
of strain, fl ~ e x (E • e), with a sign that depends on the geometry of cells (typically 
positive for prolate cells, and negative for oblate cells) [60, 68]. As was the case for the 
attraction/repulsion between cells, both oj and E change sign with p, and qualitatively 
different rotational behaviors are expected to occur for pushers vs. pullers. Nearby 
pushers on a converging course induce flow fields on each other that reorient them 
in the side-by-side configuration (Fig. 11c). If the cells are separated by a distance 
h, and are oriented with an angle 9 with respect to the distance perpendicular to 
their separation, the reorientation takes place with a rotation rate ft <~ —p9/nh 3 . In 
contrast, pullers induce the opposite flow field, which leads to a reorientation of the 
cells in the elongated direction (Fig. lid). Interestingly, for both pushers and pullers, 
the final configuration is one for which cells induce attracting flow fields on each other. 
As discussed below, similar results govern the orientation of cells near boundaries. 

The results above describe appropriately the leading-order hydrodynamic 
interactions between cells. Higher-order effects can be considered with various levels of 
modeling accuracy, two of which we address here. First, there is an active component 
to hydrodynamic interactions. The physical reason is the following: In the flow 
field created by cell A, cell B sets up its own disturbance flow (just as a solid body 
sets up a disturbance flow when located in a shear flow), which then influences the 
velocity and orientation of cell A. This is a weaker effect (decaying faster in space), 
but important for artificial swimmers with no permanent dipoles [174]. Higher-order 
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"reflections" can be considered in a similar manner. In addition, the flow induced by 
a cell is only dipolar at leading order, and includes also in general higher singularities 
(decaying faster), such as source-dipolcs and force quadrupoles. The quadrupolar 
contribution is important because it dominates the pair velocity correlations, as 
recently demonstrated experimentally for swimming E. coli [175]. 

Beyond the simple physical picture presented above, a few studies have looked 
in detail at hydrodynamic interactions between more realistic models of swimming 
cells. Interactions between two cells were studied analytically and numerically 
for two different cases, spherical squirming (swimming motion where all surface 
deformation occurs tangentially to the swimmer surface) [176, 177] and swimming 
with a single helical flagcllum [178]. Detailed experiments were also carried out 
to study hydrodynamic interactions between the protozoan Paramecium [177]. In 
all cases, the far-field physical picture described above is correct, but the details of 
hydrodynamic interactions at short range are also important, and in many cases lead 
to an instability of the side-by-side configuration for pushers, and unsteady three- 
dimensional cell trajectories. 

Beyond the dilute limit, dense cell suspensions display remarkable complex 
dynamics. Oriented suspensions of swimmers are found to be long-wavelength unstable 
[179, 180], with persistence of short-range order [180], leading to nonlinear states of jets 
and swirls [181] similar to those observed experimentally [161, 162, 163, 164, 165, 166] 
(see Ref. [182] for a two-dimensional study). Isotropic suspensions have also 
been found to be unstable [181]. As a result of the large-scale motion in dense 
suspensions, swimming cells undergo a remarkable decrease in effective diffusivity 
[183, 184, 185, 186]. 

One final important feature of cell suspensions is their rheological characteristics 
[187, 188]. In the limit where cells do not interact with each other hydrodynamically, 
the response of a cell population to an external shear can be quantified using 
Batchelor's theory for suspensions of force-free bodies [67]. For simplicity, let us 
consider spherical cells with radius a, described by force dipolcs p = pe, and 
distributed with volume fraction c « 1 in a Newtonian fluid of viscosity rj. If an 
external shear flow is applied to the suspension, with shear rate 7, the effective 
viscosity of the population, r] e s, is different from the background viscosity of the 
Newtonian fluid at leading order in c by the amount 



In Eq. (55), the coefficient 5/2 is the Einstein contribution to the viscosity [65], 
(...) denotes averages over the cell population, the directions 1 and 2 refer to the 
external flow and shear directions respectively, and t s = pja? is the typical active 
stress created by the swimming cell. For swimming E. coli, p w 0.1-1 x 10~ 18 Nm 
[189] and a w 1-10 /zm, so t s rs 10 _4 -1 Pa. Since the direction e of each cell will 
rotate as a result of the external shear, the effective viscosity of the population, 
Eq. (55), is unsteady and shear-dependent. When t s <C 777, the viscosity is 
dominated by the passive (Einstein) response of the cells, whereas when t s ^> 777 
the rhcology is expected to be governed by active stresses. In addition, anisotropy in 
the distribution of swimming cells (e.g. orientationally ordered state [187]) leads in 
general to normal stress differences with normal stress coefficients, vPi and *2, scaling 
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A detailed study of the rheological characteristics of non-dilute suspensions of cells 
was carried out for spherical squirmers [188]. Hydrodynamic interactions give rise to 
time- varying cell-cell configurations, and affect therefore the rheological properties at 
order 0(c 2 ) [67]. 

7.2. Interactions between cells and boundaries 

Just as other nearby cells influence the dynamics of a swimming microorganism, 
the presence of boundaries, and more generally confinement, significantly impacts 
cell locomotion. In addition to affecting the concentration of chemical species that 
influence the motility of microorganisms [199], boundaries modify the hydrodynamic 
stresses acting on the cells, and near-wall motility is therefore both biologically 
and physically different from bulk motility. Biological locomotion near boundaries 
includes surface-associated bacterial infections [200, 201], biofilm formation [202, 203], 
spermatozoa locomotion at the uterotubual junction [2] and surface-associated 
behavioral change [204] . In this section, we focus on the fluid mechanics of locomotion 
near walls. 

Four distinct aspects of cell locomotion are modified by the presence of nearby 
boundaries. The first is the change in the swimming speed near a wall, which was 
addressed theoretically by a number of studies [94, 205, 206, 207, 13]. Since viscous 
drag increases as a body comes closer to a boundary, it might be expected that a cell 
would slow down. However, since the propulsion method is also drag-based, a closer 
look is necessary. Indeed, for the Taylor sheet geometry, we mentioned in §5.1 that 
the presence of a nearby wall speeds up the swimmer with a prescribed waveform. 
Here consider this matter further for a swimmer that has no head, and swims using 
planar waves on a flagellum. As was shown in Eq. (35), for a given waveform of 
the flagellum, the swimming speed of the cell is an increasing function of the ratio 
between the perpendicular and parallel (£||) drag coefficients. Both coefficients 
are found to increase near a wall, but £j_ increases faster than £|| , so that the ratio 
increases, and so does the swimming speed. Physically, for fixed waveform, the 
drag-based propulsive force generated by the flagellum increases near the wall, and so 
does the resistive drag on the swimmer, but the propulsion increase is stronger, and 
therefore the swimming speed goes up. Associated with the increase in the speed, 
there is an increase of the rate of working that the swimmer has to provide in order to 
maintain the same waveform near the wall [94, 205, 206, 207, 13, 208]. If alternatively 
the swimmer is assumed to swim with constant power, then the presence of a boundary 
leads in general to a decrease of the swimming speed, leading to a decrease of the 
swimming efficiency except for some special flagellar waveforms [94, 205]. 

The second type of wall-influence on the swimming kinematics of some 
microorganisms is a change in their trajectories. This is most famously the case 
for swimming bacteria which have helical flagclla, such as E. coli, and change their 
swimming trajectory from straight to circular near a surface (Fig. 12a). In an infinite 
fluid E. coli swims, on average, in a straight line, with a swimming speed found by 
balancing flagellar propulsion by drag (see §5.2). The chiral shape of the flagella is 
important for propulsion generation, but since the propulsive force is axisymmetric 
when averaged over on period of flagella rotation, the motion occurs on a straight 
line. Near a wall, the chiral propulsion mechanism leads a breaking of the (time- 
averaged) axisymmetry of the propulsive force because new non-zero components arise 
in the motility matrix of the helix. For a helix parallel to a surface, the presence of 
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(b) 



Figure 12. Circle swimming of bacteria near boundaries. (a): Circular 
trajectories for smooth-swimming E. coli bacteria near a glass surface; 
superposition of 8 seconds of bacteria position (picture courtesy of Willow 
DiLuzio); (b): Physical mechanism for circle swimming; the rotation of the 
bacterial fiagella near the surface (blue dotted arrow) induces a net force on the 
flagella parallel to the surface but perpendicular to the fiagella axis (red arrows, 
solid); an equal and opposite force is acting on the cell body due to its counter- 
rotation, resulting in a wall-induced torque acting on the cell, and a circular 
swimming trajectory (black arrow, dashed). 



a boundary leads to a non-zero coupling between the rotation of the helix around 
its axis (Fig. 12b, blue dotted arrow), and the force on the helix in the direction 
perpendicular to the helix axis and parallel to the surface (Fig. 12b, red solid arrow). 
In other words, when the helical fiagella rotate, they create a net force on the cell 
at a right angle with respect to the motion and parallel to the surface. There is 
an exact and opposite force acting on the cell body, which rotates in the opposite 
direction as the flagella, and the net effect is a wall-induced torque (Fig. 12b, red 
solid arrows). If the cell were to continue swimming in a straight line, it would have 
to apply a net torque on the surrounding fluid. Since a swimming bacterium is in 
fact torque-free, the cell cannot swim straight but instead rotates at a rate such that 
the viscous torque from that rotation exactly balances the wall-induced torque, and 
therefore swims along circles on the surface (Fig. 12b, black dashed arrow). For a cell 
using a left-handed helix for propulsion, such as E. coli, this effect leads to swimming- 
to-the-right [209, 210, 211, 212, 213, 214]. 

The third influence of boundaries on swimming cells is analogous to the attraction 
and re-orientation induced by cell-cell interactions and discussed in §7.1. Consider a 
single cell moving near a solid wall. As the cell is swimming, it sets up a dipolar flow 
field, but this flow field in general does not satisfy the no-slip boundary condition on 
the wall. The flow in the presence of the wall is then the linear superposition of the 
infinite-fluid flow, plus any image flow field created on the other side of the surface, 
necessary to enforce the boundary condition on the wall exactly. Image flow fields for 
the fundamental singularities of Stokes flows were pioneered by Blake [215, 216], and a 
textbook treatment is offered in Refs. [60, 63]. Physically, the approach is similar to the 
method of images in electrostatics, but with the complication that in fluid mechanics 
the no-slip boundary conditions is vectorial, whereas the constant-potential boundary 
conditions in electrostatics are scalar. Because of the image flow field, a cell described 
by a dipole strength p, located at a distance h from the surface, and pointing at an 
angle 9 from the surface direction (see notation on Fig. 13a) is subject to the gradients 
of the image flow field, and as a result rotates with speed fl ~ —pd/nh 3 in the direction 
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Figure 13. Wall-induced rotation of swimming cells. A swimming cell is located 
at a distance h from a solid surface, and at an angle 8 with respect to the direction 
of the surface, (a): Pushers are reoriented hydrodynamically in the direction 
parallel to the surface (equilibrium, 8 = 0); (b): Pullers are reoriented in the 
direction perpendicular to the surface (equilibrium, 8 = ±tt/2). 



parallel to the surface and perpendicular to the cell body [189]. The l/h 3 scaling 
originates from the leading order vorticity of the image flow field, which is also dipolar. 
The rotation occurs as if the cell is interacting hydrodynamically with a mirror-image 
cell located on the other side of the surface, and the rotation rate is therefore analogous 
to that quantified by Eq. (54) §. If the cell is a pusher, the wall-induced rotation rate, 
£1 = d8/dt, tends to align the swimming cell in the direction parallel to the surface 
(9 = 0, Fig. 13a). As a result of this parallel configuration, the cells will swim in 
a side-by-side configuration with their image cell, and are therefore attracted to the 
surface with an attractive speed scaling as u± ~ p/ijh 2 (see Eq. 53). This physical 
picture explains the accumulation of swimming cells near surfaces observed in many 
biological experiments [217, 218, 219, 208, 220, 221, 183, 189]. In contrast, cells which 
are pushers are rotated in the opposite direction, and their stable configuration is 
instead at a right angle with respect to the surface (9 = ±ir/2, Fig. 13b). In a confined 
environment, these cells are therefore always swimming toward one surface, a result 
which also leads to their accumulation. The physical origin is however different, and 
instead of being attracted by an image cell as pushers are, pullers simply swim into 
the wall. 

The fourth hydrodynamic effect of boundaries, less studied, is a potential 
reduction of cell-cell hydrodynamic interactions near solid surfaces. Indeed, in many 
cases, a flow singularity at a distance h from a solid surface is canceled out in the far 
field by its image on the other side of the surface, and the overall flow decays faster 
in space than the original singularity (at distances further away than h from it). For 
example, a stokeslet near a wall has a dipolar behavior in the far field if the stokeslet is 
parallel to the surface, and a quadrupolar decay if it is perpendicular to it [215, 216]. 

§ The effect of a flat boundary is mathematically equivalent to the presence of a mirror-image cell 
instantaneously located on the other side of the surface if the surface is a no-shear interface (e.g. a 
free surface). If instead it is a no-slip surface, the analogy is not quite exact mathematically, but it 
remains qualitatively correct. 
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Swimming cells typically behave as force-dipoles in an infinite fluid, but when near a 
solid surface, configurations exist where the total flow field (dipole + images) ends up 
decaying as 1/r 3 or 1/r 4 , and hydrodynamic interactions with other cells are weaker 
as a result. In that case, the distance h between the swimming cell and the wall acts 
therefore as an effective cut-off for cell-cell interactions. 

7. 3. Interactions between flagella 

7.3.1. Eukaryotic flagella and cilia. We now decrease our length scales, and 
discuss hydrodynamic interactions between cellular appendages. The first historical 
evidence of important hydrodynamic effects between nearby eukaryotic flagella were 
observations by James Gray of in-phase beating of spermatozoa flagella, later 
reproduced by Taylor [29]. Similar observations have been made since then, such as 
the synchronous beating of Fishfly spermatozoa [170]. Why would organisms phase- 
lock in such a manner? Exploiting the swimming sheet model, Taylor showed that two 
swimmers with the same prescribed waveform (Fig. 14a) dissipate the least amount of 
mechanical energy for swimming when the waves are exactly in phase (Fig. 14b), and 
the dissipation goes up monotonically with phase difference. Using high-amplitude 
numerical simulations, these results were later revisited [222]. Two sheets with same 
waveform but different phases are seen to swim with different velocity, and their phase 
difference 4> evolves until they are either perfectly in phase (<fi — 0) or perfectly out of 
phase (4> — 7r), with both configurations being stable [222]. 

Two nearby swimmers therefore display phase-locking, but the locked state can 
be that of maximum dissipation, in contrast with what is observed experimentally. 
For real eukaryotes, the flagellar waveform is however not fixed but is found as a 
solution to a mechanical problem, as explained in §6: The shape arises from a balance 
between internal forcing inside the flagella (molecular motors), passive elastic (possibly 
viscoelastic) resistance from the axoneme, and viscous resistance from the outside 
fluid. When there is a second swimmer located nearby, the fluid forces are modified, 
and so is the mechanical force balance determining the flagella shape. Through this 
neighbor-induced change in the waveform, cells arc expected to be able to perfectly 
phase lock, as is observed in the synchronization of nonlinear oscillators [223] . 

Most of the work related to hydrodynamics-induced waveform change and 
phase locking was concerned with cilia, which are short eukaryotic flagella (a few 
microns long), found in densely packed arrays and involved in many important 
biological functions, including sensing, fluid transport, locomotion, and development 
[224, 13, 3, 58, 21, 225, 226]. An individual cilium deforms in a non-reciprocal fashion, 
with a typical high-friction power strokes (sequence 1 to 3 in Fig. 14c) and a low friction 
recovery stroke (sequence 4 to 8 in Fig. 14c). The internal actuation of each cilium 
is independent from that of its neighbors, and they do not communicate with each 
other except through the fluid. However, when they are closely-packed on surfaces, 
cilia arrays display collective behavior, termed metachronal waves. The deformation 
of each cilium is locked in phase with that of its neighbor, with a constant (small) 
phase difference, leading to propagating waves of deformation (Fig. 14d) [58]. 

The physical origin of coordinated beating is the central question of ciliary 
dynamics, which a number of theoretical studies have attempted to answer [108, 109, 
110, 111, 227, 228, 229, 230]. Two different approaches have been proposed. In the first 
one, the mechanics of each cilium is modeled in the most accurate way, and numerical 
simulations are used to compute the collective beating [108, 109, 110, 111]. The 
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Figure 14. Hydrodynamic interactions between eukaryotic flagella and cilia, 
(a): Nearby swimming cells with arbitrary phase difference; (b) Phase- locked 
nearby swimming cells; (c) The asymmetric beat pattern of a cilium (sequence 1 
to 10); (d) Schematic representation of the metachronal waves of collective beating 
displayed by arrays of cilia. 



crucial ingredient in that approach is to correctly model the internal load-dependent 
force generation in the axoneme (without load dependence, or feedback, there is no 
phase- locking as can be expected [230]). With that approach, it is found that two 
cilia starting randomly end up beating in perfect synchrony within two beating cycles 
[110]. If instead there is a large number of cilia, waves arise naturally as a result of 
hydrodynamic interactions [110]. Subsequent work showed that as the waves develop 
from arbitrary initial conditions the rate of work done by the cilia as they are beating 
is decreasing [111]. Physically, because of viscous drag, it is energetically advantageous 
for one cilium to beat in the presence of a neighboring cilium with similar phase. 

The second approach considers simplified models for the dynamics of the cilia, 
providing analytical insight into the necessary conditions for phase-locking [227, 228, 
229, 230]. A first study considered a regular lattice of cilia, where their direction 
(the direction of the beating plane) is assumed to obey a balance between rotational 
Brownian motion and rotation induced by the flow created by all other cilia. For 
small enough temperature, a transition is observed between a state where no net flow 
occurs on average, to a state where all cilia point in the same direction and a net 
flow is created [229]. Further insight is provided by considering a simplified load- 
dependent internal molecular engine. In that case, metachronal waves arise only if a 
constant phase shift is assumed to exist between each cilium and its neighbor [229]. 
Motivated by nodal flows in development [226] , a second study considered cilia whose 
tips perform three-dimensional trajectories over a surface. Each cilium is modeled by 
a sphere subject to an active load-dependent force, and interacting hydrodynamically 
with a second cilium. Depending on the relative position and orientation of the two 
cilia models, in-phase (<f) — 0) or out-of-phase locking (<f) — n) arise from random 
initial conditions [228]. A similar model with two sphere- like cilia rotating due to 
an applied torque near a wall was recently proposed. In that case, in-phase locking 
is obtained provided that the circular trajectory of each cilium is allowed to vary in 
response to hydrodynamic interactions [230] . 



7.3.2. Bacterial flagella. Hydrodynamic interactions between flagella also play a 
pivotal role for bacterial locomotion. In that case, the phenomenon of interest 
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Figure 15. Bundling of bacterial flagella. During swimming, the bacterial flagella 
are gathered in a tight bundle behind the cell as it moves through the fluid (a 
and d). During a tumbling event, the flagella come out the bundle (b), resulting 
in a random reorientation of the cell before the next swimming event. At the 
conclusion of the tumbling event, hydrodynamic interactions lead to the relative 
attraction of the flagella (c), and their synchronization to form a perfect bundle 
(d). 



is flagellar bundling [231, 232]. Wild- type swimming bacteria, such as E. coli, 
typically display "run-and-tumble" behavior during their locomotion. During runs, 
the bacterium swims along a roughly straight path, and its flagellar filaments are 
bundled together tightly behind the cell (Fig. 15a and d). Near the end of a run, one 
or or more motors reverses, and the corresponding filaments unwind from the bundle. 
Viscous stresses lead to polymorphic transitions which turn the cell, ultimately leading 
to a random reorientation once the motors reverse again and the full bundle forms [22] 
(Fig. 15b). The process by which all the flagella, which are randomly distributed along 
the cell's surface, come together (Fig. 15c) to form a perfect aggregate at the end of a 
tumbling event (Fig. 15d) is termed bundling. Since the flagella do not communicate 
except through the fluid, it has long be postulated that bundling occurs only through 
hydrodynamic interactions. 

During bundling, two different physical mechanisms are involved: (1) the 
attraction between the rotating flagella, and (2) their subsequent phase-locking inside 
the bundle. Concerning the attraction between the rotating flagella, two different 
scenarios have been proposed. Both mechanisms play a role. One is purely kinematic, 
and relies on the simple observation that as the cell starts swimming, the drag on all the 
flagella naturally sweeps them behind the cell [231]. Under this scenario, the flagella 
are not so much attracted by each other, but are simply passively dragged behind 
the cell body. The wrapping of the filaments in a helical shape is then achieved 
in a passive manner by the body rotation, which is necessary in order to achieve 
overall torque balance (see §5.2) [231, 233]. The second attraction scenario relies on 
hydrodynamic interactions between the flagella, in which each helix induces a flow 
which causes the other to bend and twist about its neighbors. The geometry of the 
helices is critical for this scenario. The role of this geometry has been studied with 
macro-scale experiments using flexible helices in air [234]. In these experiments, the 
bundling was forced by means of guides at the distal ends of the helices. It was found 
that left-handed helices rotating counter-clockwise as viewed from the distal end (the 
same handedness and rotation sense as a bacterium on a run) can maintain a steadily 
rotating bundle, without jams, when the helices are wrapped around each other in a 
right-handed sense. Later, these experiments were extended with polymeric helices in 
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highly viscous silicon oil [235, 236]. There it was shown that induced flow naturally 
causes tight bundling of pairs of left-handed helices driven counter-clockwise. The 
mirror image of this arrangement but no other combination also led to bundling [235]. 
The time for bundling in the experiment was observed to scale with the rotation 
speed of the helix, rather than the relaxation time of flexible helices in a viscous fluid. 
Bundling has been also been studied computationally using a regularized version of 
slender-body theory [237]. For three left-handed helical flagella driven at constant 
torque and rotating counter-clockwise, hydrodynamic interactions lead to a bundling 
of the filaments, although not as tight a bundle as observed experimentally [237]. 

The issue of synchronization between helices driven with constant torques was 
addressed numerically [238, 239]. Perhaps surprisingly, two nearby rotating helices are 
seen to not synchronize if the helices are perfectly rigid [238], and some mechanical 
flexibility is necessary for flow-based synchronization to occur [239]. In bacteria, this 
flexibility is likely to be provided by the flagellar hook which connects the rotary motor 
to the helical flagellum [240, 241]. The contact forces that arise between different 
helices as the bundle forms likely also play a role in synchronization. 

8. Swimming in complex fluids 

Biological fluids are often laden with polymers and therefore have non-Newtonian 
rhcology. For example, mucus is found at the cervix and throughout the rest of the 
mammalian female reproductive tract [2] . The nature of the cervical mucus influences 
the likelihood of fertilization in humans: There is a good correlation between the 
hydration of the cervical mucus and the incidence of pregnancy [242]. Although 
cervical mucus is a complex, heterogeneous gel, some simple trends have been observed. 
There is an inverse relation between the degree of viscoelasticity of the mucus and 
the ability of sperm to penetrate the mucus [243, 244, 245, 246]. Human sperm 
flagella in cervical mucus have a higher beat frequency, smaller amplitude, and shorter 
wavelength than when in semen [247] . The swimming speed is the same in both media, 
but the sperm swim along straighter paths in cervical mucus [247] . In addition to the 
mucus at the cervix and in the uterine cavity, sperm must also penetrate the matix 
coating the ovum, known as the cumulus oophorus, which is a viscoelastic actin-based 
gel [248]. Spermatozoa can undergo an internal change known as hyperactivation in 
which the beat pattern changes from a symmetric to asymmetric form; hyperactivated 
cells have better motility in viscoelastic media [249, 250, 251, 252, 253]. 

Rheological measurements show that cervical mucus is highly viscous, with a 
viscosity of 100 Pas [254]. The viscoelasticity is characterized by a relaxation time of 
r = 1-10 s [255, 256, 257, 245] and an elastic modulus of 0.1-10Pas [257, 245]. These 
properties depend strongly on the phase of the menstrual cycle as well as hydration. 
Since the typical beat frequency of a flagellum is uo = 20-50 Hz [13], the Deborah 
number lut is easily larger than one, indicating that viscoelastic effects are important. 

There are several other important examples of swimming or transport in complex 
fluids. The cilia that line the human upper airway lie in a thin layer of of Newtonian 
liquid, which in turn is coated with a high- viscosity layer of mucus [58]. Again, 
since the Deborah number for this system is large [258], elastic effects are important. 
Bacteria also encounter viscoelastic fluids: The ulcer-causing bacterium Helicobacter 
pylori swims through mucus lining the stomach [259], and spirochetes move through 
the connective tissue of the host during infection [260] . 

In our discussion of swimming in a complex fluid, we use continuum mechanics 
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to model the fluid. This simplification is necessary for making progress, but it is 
important to realize that the size scale of the microstructure of the fluid can be 
comparable to the size of the swimmer, and therefore a different approach may 
be necessary to accurately capture the interactions between the swimmer and the 
polymers. Complex fluids display a vast array of non-Newtonian effects, such as stress 
relaxation, normal stress differences, and shear-rate dependent viscosity [191, 192, 
193, 194, 195, 196, 197, 198]. Our approach is to illustrate some of the distinctive 
properties of swimmers in complex fluids by focusing on one class of models, fluids 
with fading memory. These models apply to polymer solutions. 

When a polymer solution is subject to shear, the polymers stretch out. The 
resulting loss of entropy of the polymers leads to an effective elastic force tending to 
recover the initial configuration of the polymers. The balance of the entropic force and 
the viscous drag on the polymer sets the time scale A over which the fluid has memory. 
The appearance of this new time scale gives polymeric liquids a completely different 
character relative to viscous Newtonian fluids. First, the property of kinematic 
reversibility is lost, even if the Reynolds number is vanishingly small, since growth or 
decay of stress can lag the change in shape of a swimmer by the time scale A. Second, 
the new time scale implies that the constitutive relation for a polymer solution depends 
on the rate of change of stress with time, which automatically leads to nonlinear 
terms. This fact may be seen by invoking material frame indifference, a fundamental 
assumption of continuum mechanics [190, 192], or by deriving the continuum theory 
directly from a microscopic theory of the polymers in solution [196, 198]. Many of 
the properties of Stokes flow that we have invoked in our study of swimming have 
relied on the linearity of Stokes flow. The nonlincarity of the equations of motion for 
a polymeric liquid generally implies that the scallop theorem does not hold. The loss 
of linearity spoils superposition, and therefore we may not use the singular solutions of 
§3.3 to construct a slender body theory. Hydrodynamic interactions between distant 
cells will have a different character than in Stokes flow, since the far field form of the 
velocity field due is changed. These features of polymeric liquids make the study of 
swimming in complex fluids a daunting challenge, but also an area of opportunity. 
Except for a few early works [261, 262, 263], the area is largely unexplored. 

To examine these issues, we review the calculation [258] of the swimming velocity 
of the Taylor sheet of §5.1 for a fluid described by the Oldroyd-B model [192]. In 
this model, the deviatoric stress r = pi + er is the sum of a contribution from the 
polymer solute t p and a contribution from the Newtonian solvent r s = 2?y s E, where 
E = ^(Vu + Vu T ). The polymer contribution relaxes to the viscous stress over the 
time scale Ai, 

-r p + Ai t p = 27? P E, (57) 

where 

t= — ^ + u • Vr - t • Vu - (Vu) T • r (58) 
at 

is the upper-convected derivative of r. The upper-convected derivative of a tensor r 
is the expression in general coordinates of the rate of change of the tensor calculated 
in a frame that translates and deforms with the local fluid velocity [190, 192]. Using 
t = t p + t s and eliminating r p from Eq. (57) yields the Oldroyd-B model, 

t + Xi t= 2i](E + A 2 E) (59) 



The hydrodynamics of swimming microorganisms 



42 



where 77 = f] p + rj s and A 2 = ffc^i/ 7 /' < The equation of motion is —Vp + V • r = 
with V • u = 0. Note that the explicit presence of the time derivatives as well as the 
nonlinear terms in Eqn. (59) spoils kinematic reversibility and violates the assumptions 
of the scallop theorem. 

The solution to the swimming problem for a sheet with prescribed traveling wave 
(20) proceeds just as in the Newtonian case. The velocity field, stress, and boundary 
conditions are expanded to second order in amplitude b, and the governing equations 
are solved subject to the boundary conditions (21) and (22). To first order, the flow 
field is identical to the flow field of the Stokes problem. However, due to the relaxation 
time A, the first-order component of the stress field r has a lag relative to the Stokes 
stress field. Through the nonlinear terms of Eq. (59), this stress field drives a second- 
order flow that leads to the swimming speed 

itti 1 7.12 1 +u 2 \ 2 1 rj s /v 

W = r kb i+^a; • (60) 

Since rj s < 77, the swimmer in an Oldroyd-B fluid is slower than in a Stokes fluid. Note 
that unlike the Stokes, case, the swimming speed (60) depends on material parameters 
such as Ai , rj s , and 77. The Oldroyd-B constitutive equation is inapplicable to flows with 
large extension rates. Nevertheless, the result (60) continues to hold for more accurate 
fading memory models such as FENE-P, the Johson-Segalman-Oldroyd model, and the 
Giesekus model [258]. The expression for the speed, Eq. (60), also continues to hold 
for a cylindrical filament with a traveling wave, in the limit in which the radius of the 
cylinder is much smaller than the lateral displacement of the cylinder [148, 264] . 

The nonlinear dependence of the swimming speed on the frequency signals the 
breakdown of the scallop theorem. To see why, first consider the Newtonian limit 
Ai — > 0, in which all elastic effects vanish and the speed is proportional to u>. Now 
consider a reciprocal motion situation in which a traveling wave of wave number k and 
frequency uj travels rightward on the sheet for a period (27r/w)/3, and then leftward 
for a period 2(2it/uj)/3 [264]. The sheet will have zero net displacement after this 
process. Now consider this waveform for a sheet in an Oldroyd-B fluid. Suppose the 
sheet has periodically been executing these motions long enough that transients from 
startup from rest have died away. Since the speed now depends nonlinear ly on u), the 
translation in each segment of the motion is different: The sheet moves faster during 
the leftward motion with the smaller frequency. Thus, there is a net displacement. 
This argument disregards the memory of the fluid, since the stress during each stroke 
is effected by the previous stroke. Accounting for these memory effects leads to a 
slightly smaller net displacement, but it still has the same qualitative dependence on 
u)X\. 

The properties of the medium can affect not only the speed of a swimmer with a 
prescibed stroke, but also the shape of a beating flagella. For example, as the viscosity 
of the solution increases, the wave form of a human sperm flagellum flattens along most 
of its length, with most of the deflection taking place at the distal end [265]. Thus it 
natural to model the dependence of the shape on material parameters such as viscosity 
and relaxation time. As we learned in §6, the small- amplitude shape is determined by 
an equation which is first-order in amplitude. Thus, we only need the fluid force to 
first order. For the Oldroyd-B fluid, the calculations for the filament geometry shows 
that this fluid force has the form of resistive force theory, with complex frequency- 
dependent effective viscosity [87] . Denoting Fourier components with a tilde, the shape 
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Figure 16. Artificial swimmers with discrete degrees of freedom: (a) Purcell's 
three-link swimmer (red arrows indicate change in conformation of the swimmer 
shape); (b) Macro-scale experimental realization of Purcell's swimmer (Photo 
courtesy of Brian Chan and Peko Hosoi); (c) Three-sphere swimmer; (d) Micro- 
scale realization of a three-sphere pump [267] ; (e) Two-sphere swimmer with shape 
change; (f) Hydrodynamic interactions between two reciprocal swimmers. Picture 
in Fig. 16d adapted from Ref. [267] courtesy of Marco Cosentino Lagomarsino. 



equation for an Oldroyd-B fluid is 



i£-L 7 ^h df 

: h = - A ^z~ a ^ ( 61 ) 



1 + iAia; dx dx ' 

where / is the sliding force density. By prescribing the force density, imposing the 
boundary conditions, and using linear superposition (valid to first order), we may 
calculate the shape of the filament as a function of material parameters. The result 
is that this model gives patterns that look qualitatively similar to the experimental 
observations [87] . Since the shape of the filament determines the swimming velocity, 
the dependence of shape on relaxation time gives an addition correction to the 
swimming velocity of the prescribed sheet. 

A natural next step is to extend the study of swimming to large amplitude 
deflections. Recent numerical work on large-amplitude peristaltic pumping of an 
Oldroyd-B fluid two dimensions has shown similar qualitative differences with the 
Newtonian case, such as reduced pump rate, dependence on material parameters, and 
loss of kinematic reversibility [266]. The application of numerical methods such as 
these to swimmers will be crucial for understanding finite-size effects such as the role 
of the shape of the cell body or the end of the flagellum, since large normal stresses 
can develop at such regions of high curvature. 

9. Artificial swimmers and optimization 



In this final section, we look beyond the biological realm and survey the bioengineering 
advances in design and optimization of artificial swimmers and bio-inspired systems. 
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(a) (b) (d) 

Figure 17. Artificial swimmers with continuous deformation: (a) Synthetic 
flagella made of actuated micron-scale magnetic filaments [279]; (b) Macro- 
scale swimming by boundary actuation of an elastic filament [280]; (c) Toroidal 
swimming; (d) Swimming by surface treadmilling. Picture in Fig. 17a adapted 
from Ref. [281] courtesy of Remy Dreyfus and by permission from Macmillan 
Publishers Ltd, copyright 2005. 



9.1. Designing artificial swimmers 

A number of designs have been proposed and tested for artificial swimming devices 
at low Reynolds number. They can be sorted in three different categories. The first 
category includes all swimmers which deform their shape with only a finite number of 
degrees of freedom, actuated in a non-reciprocal fashion (and therefore, at least two 
degrees of freedom are needed). The original example is Purcell's three link swimmer, 
which posseses two hinges varying in time with phase differences (Fig. 16a and b) 
[14, 70, 268, 75]. A second related design is that of a three-sphere swimmer (Fig. 16c): 
The distances separating the spheres vary in time and with phase differences, leading to 
locomotion [76, 269, 172, 270, 271, 272]. The extension to N spheres was also proposed 
theoretically [273]. The three-sphere design was implemented experimentally using 
optical tweezers, and the out-of-phase motion of the three colloidal beads leads to 
fluid pumping (Fig. 16d) [267]. A third design only needs two spheres, but is coupled 
with a change in shape of one of the spheres, leading to non-reciprocal kinematics 
(Fig. 16e). As the sphere increases in size, it serves as an anchor against which 
the swimmer can push, and rectify a time-periodic change in size [274]. A related 
idea using rotating two-sphere swimmers, was also proposed [275]. Finally, because 
there is no "many-scallop" theorem [276], two reciprocal non-swimmers can exploit 
hydrodynamic interactions to swim collectively (Fig. 16f). The resulting collective 
swimming speeds depend on the separation distance between the swimmers as well as 
their relative position and orientation [277, 278]. 

The second category of artificial swimmers includes bodies deforming in a 
continuous fashion. In an experimental breakthrough, a team has been able to exploit 
the response of paramagnetic filaments to external magnetic fields to construct the 
first artificial micron- long swimmer (Fig. 17a) [279], thereby motivating a number of 
theoretical studies [282, 283, 284, 285, 286]. In this work, a AC magnetic field is 
applied in the direction perpendicular to the filaments, and the presence of a body 
(here, a red blood cell), breaks the right- left symmetry in the continuous response 
of the filament. The result is a wave of deformation traveling from the tip of the 
filament to the body it is attached to, and locomotion flagella-first [279]. A second 
type of swimmer with continuous deformation exploits the boundary actuation of 
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Figure 18. Locomotion powered by chemical reactions [281]. Polystyrene 
colloidal spheres in a solution of water and hydrogen peroxide ([H2O2]) with 
increasing concentration. Top: Control experiment (inert spheres); Bottom: 
Experiment where the spheres are half-coated by platinum, a catalyst for the 
reduction of [H2O2] into oxygen and water, which leads to directed swimming 
of the spheres (coupled to Brownian rotational diffusion). Picture reproduced 
from Ref. [281] courtesy of Ramin Golestanian by permission from the American 
Physical Society, copyright 2007. 



elastic filaments to generate propulsive forces (see §6), and locomotion (Fig. 17b). In a 
macro-scale experiment, the forces arising from such elastic swimming were measured 
[130], and the actuation method was implemented successfully to obtain swimming 
[131, 280]. A third type of swimmers with continuous deformation are those with 
a special mode of surface deformation called "surface treadmilling" , for which the 
shape of the swimmer is fixed and where its body undergoes a continuous tangential 
displacement along its surface (Fig. 17c, d). In his original article, Purcell proposed 
such design in the shape of a torus undergoing continuous surface rotation (Fig. 17c) 
[14], an idea recently analyzed in detail [287, 288, 289]. If the body is slender and 
displays directed tangential displacement all along its surface, so that a material source 
is present on one side of the body, and a sink on the other side (Fig. 17d), the resulting 
locomotion can occur with arbitrarily high efficiency [290]. 

The final category of artificial swimmers use chemical reactions to power 
locomotion, a case investigated both experimentally [291, 292, 293, 294, 295, 296, 
281, 267, 297] and theoretically [298, 299, 300]. The prototypical setup is a body 
composed of two materials, one that is inert, and one that is a catalyst or reactant 
for a chemical reaction. The presence of a chemical reaction leads to an imbalance in 
concentration for some of the reactants and/or products to the reaction, which leads 
to an imbalance of osmotic pressures on the body, and results in directed propulsion. 
An example of the trajectory of such chemical swimmer is illustrated in Fig. 18. 



9.2. Exploiting low-Re locomotion 

The physical mechanisms of propulsion at low Reynolds number can be exploited 
in a number of ways beyond biological or synthetic locomotion. In a pioneering 
experiment, suspensions of swimming bacteria were seen to provide an effective high- 
temperature thermal bath for suspended inert particles [162]. Similar results were 
obtained for diffusive mixing of two fluids containing swimming cells [301], and for the 
motion of colloidal particles above surfaces covered with attached bacteria (Fig. 19) 
[302] . A second utilization of swimming is cargo-towing, with potential applications 
in biomedical devices. Experiments were performed with solid bodies covered with 
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Figure 19. Enhanced diffusion near a carpet of swimming cells [302]. Trajectory 
of fluorescent particles over a 10 s interval when located above a flat surface 
coated with attached bacteria (two-dimensional view from above the surface). 
(A): The particles are located 80 fim above the surface, and move by Brownian 
motion. (B): When the particles are located 3 (im above the surface, their effective 
diffusion is significantly larger than Brownian motion [302], Picture reproduced 
from Ref. [302] courtesy of Kenny Breuer and by permission from the Biophysical 
Society, copyright 2004. 



attached bacteria, with successful translational and rotational towing [302] . A related 
method was also devised to transport micron-size objects at the single-cell level 
[303], and a theoretical framework for towing by swimming now exists [304]. A final 
application of flagella motion is pumping. Indeed, swimming and pumping are dual 
problems, and in general a tethered swimmer acts as a pump [305]. One approach, 
studied numerically, uses elastic filaments attached to a solid surface and actuated 
by external time- varying torques [306] . Experiments were also conducted on surfaces 
covered with filaments made of self-oscillating gels and displaying wave-like motion 
[307]. 

9. 3. Optimization 

Motivated by the optimal tuning of synthetic micro-swimmers, as well as possible 
insight into evolutionary processes, significant work as been devoted to the 
optimization of locomotion at low Reynolds number. Since time can be scaled out of all 
low Reynolds number swimming problems in Newtonian fluids, optimizing swimming 
speeds is not, in general, a well-posed mathematical problem, and an additional form of 
normalization is required. The traditional approach is to define a swimming efficiency 
as the ratio of the useful to the total rate of working of the swimmer against the viscous 
fluid (see §5.2). The optimization problem becomes a maximization problem for the 
efficiency [11], and is equivalent to finding the swimmer with the largest swimming 
speed for a given amount of mechanical energy available. 

Two types of optimization questions can be asked, related to either the waveform 
or the overall geometry of the swimmer. The first problem was formulated by Lighthill 
as the following: If a flagellum is being distorted as a planar traveling wave, what 
would be its optimal waveform [11]? Assuming an infinitely long flagellum, and within 
the framework of the local drag theory, the optimized flagellum is a sawtooth wave, 
where each branch of the sawtooth makes an angle ±9 with the horizontal, with 
9 = arccos [(1 + (£|| /£±) 1/2 )~ 1/2 ] [11]. Since ~ 1/2, the optimal angle for the 

sawtooth is 9 ~ 40°, leading to an optimal value for the product of the wave amplitude, 
6, and the wavenumber, k, as bk w 1.3. The optimal value for 9 is insensitive to the 
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exact value of the ratio 6/^i,oron the possible presence of a cell body [1 1] . Numerical 
simulations for sawtooth patterns revealed the optimal criterion of Lighthill to remain 
valid for finite-size swimmers [308]. In addition, flagella with the optimal sinusoidal 
deformation have very similar geometrical characteristics: bk w 1.26 using the local 
drag theory [309], bk s=s 1 using slender body theory [310]. Similar results are obtained 
by numerical optimization of the flagellar waveform [311]. 

The existence of an optimal angle between the direction of swimming and the 
direction of the flagellum can be understood on the basis of the physical picture we 
introduced for drag-based propulsion (Fig. 2). Indeed, we see from Eq. (13) that 
the propulsive force arising from the local motion of the flagellum scales with its 
orientation as / pr0 p ~ sin 9 cos 9. If the work done by the flagellum was not a function 
of its orientation, the optimal angle would therefore be the one for which sin 26 <~ 1, so 
9 ~ 45°. Because £_l > some energy can be saved by decreasing the angle, thereby 
promoting tangential over perpendicular motion, which explains why the optimal angle 
is close to, but below, 45°. In the case of planar waves, the optimal flagellum waveform 
is therefore kinked (sawtooth), and the slope angle alternates between +9 and —9. In 
the case of a three-dimensional helical filament, a constant angle between the local 
tangent along the flagellum and the swimming direction can be accommodated while 
keeping a smooth waveform if (and only if) the flagellum takes the shape of a helix. 
In that case, the optimal angle was determined numerically using slender body theory 
to be 9 = 45° [312] (see also Refs. [64, 313]). 

A second design feature which can be optimized is the overall geometry of the 
swimming cell, in particular the ratio of the length L of its flagellum and the radius a 
of its body, which is assumed to be spherical for simplicity. In the case of swimmers 
exploiting planar waves, this problem can only be studied if hydrodynamic interactions 
between the body and the flagellum are properly taken into account, as otherwise 
the presence of a body always decreases the swimming efficiency (see Eq. 38). This 
problem was analyzed in detail using slender body theory for sinusoidal waveforms, 
and the optimal flagellum length to body size ratio was found to be L/a w 25 [310]. 
Subsequent work determining the optimal flagellar waveforms showed that the optimal 
ratio is close to L/a w 24 [311]. In the case of helical waves, an optimal body size 
can be determined even in the absence of hydrodynamic interactions. Indeed, we saw 
in §5.2 that cells with helical flagella need a body in order to be able to move. A 
small body prevents the swimmer from moving and is not efficient, whereas a large 
body presents too large of a drag to the overall cell, and is not efficient either: An 
optimal size must therefore exist. When hydrodynamic interactions are also taken 
into account, the optimal ratio has been determined numerically to be L/a w 10 using 
slender-body theory [312] and ~ 12 using three-dimensional boundary elements [314] 
(see also Ref. [64]). Other biologically- relevant optimization problems that have been 
addressed include the optimization of feeding current for tethered cells [315], and the 
influence of the body shape on the optimal flagellar waveform [311] and overall cell 
geometry [314]. 

Finally, work has been devoted to the optimization of simple (artificial) swimmers. 
For bodies with discrete degrees of freedom, the actuation of Purcell's swimmer 
(Fig. 16a) was optimized numerically [268, 75]. More formal mathematical work 
was devoted to the optimization of the three-sphere swimmer (Fig. 16c) using a 
control-theory framework [316]. For swimmers with continuous degrees of freedom, the 
geometrical approach to swimming at low Reynolds number was used to derive optimal 
low-amplitude swimming by surface deformation of spheres and cylinders [71, 72, 73], 
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and was later extended to large-amplitude deformation of two-dimensional bodies 
[317]. 

10. Outlook 

In this paper, our aim has been to illustrate the theoretical framework necessary to 
understand biological and synthetic locomotion at low Reynolds number. Since new 
experimental methods and setups promise to reveal the mechanisms for biological 
locomotion with ever more quantitative detail, future work in the field is likely to 
be significant. We believe in particular that there is great opportunity for theorists, 
since — as we have emphasized throughout the review — simple calculations are usually 
sufficient to gain fundamental insight into the mechanisms of locomotion. Detailed 
numerical computations will also play a crucial role, since many of the issues we 
discussed, such as hydrodynamic interactions in dense suspensions of swimmers and 
locomotion in complex fluids, involve nonlinear processes. A further challenge will be 
to integrate the understanding of basic mechanisms across multiple scales, from the 
levels of molecular motors to individual cells to large populations of cells. Finally, 
many of the implications of the ideas sketched in this review have yet to be realized 
in important areas such as the ecology of marine bacteria, the formation of bacterial 
biofilms, and the mechanics of reproduction. 
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